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Abstract 

We extend the analysis of elastic pion-nucleon scattering up to 0(p 4 ) level using the extended-on-mass-shell subtrac- 
tion scheme within the framework of covariant baryon chiral perturbation theory. Numerical fits to partial wave phase 
shift data up to -\fs = 1.13 Ge V are performed to pin down the free low energy constants . A good description of the 
existing phase shift data is achieved. We find a good convergence for the chiral series at <9(p 4 ), considerably improved 
with respect to the (9(p 3 )-level analyses found in previous literature. Also, the leading order contribution from explicit 
A(1232) resonance and partially-included A(1232) loop contribution are included to describe the phase shift data up 
to Vi = 1.20 GeV. As phenomenological applications, we investigate chiral corrections to the Goldberger-Treiman 
relation and find that it converges rapidly, and the 0(j> 3 ) correction is found to be very small: 0.2%. We also get 
a reasonable prediction of the pion-nucleon sigma term cr„ N up to 0(p 4 ) by performing fits including both the pion- 
nucleon partial wave phase shift data and the lattice QCD data. We report that <t„n = 52 + 7 MeV from the fit without 
A(1232), and cr nN =45 + 6 MeV from the fit with explicit A(1232). 

Keywords: pion-nucleon scattering, chiral perturbation theory, partial wave analysis, Goldberger-Treiman relation, 
pion-nucleon sigma term 



1. Introduction 

Pion-nucleon scattering is an important process for the understanding of chiral QCD dynamics and the interpreta- 
tion of some prominent phenomenology of strong interactionsjin. Many efforts have been made to study it. However, 
unlike the successfulness of chiral perturbation theory C^PT) 12j,[3[] in the pure meson sector, a chiral expansion in the 
pion-nucleon scattering amplitude suffers from the power counting breaking (PCB) problem in the traditional subtrac- 
tion MS - 1 schemeMj] . Many proposals have been made to remedy this problem, e.g., heavy baryon (HB) chiral 
perturbation theory JH, infrared regularization (IR) scheme @], extended on mass shell (EOMS) scheme cLIg], etc.. 

As a successful non-relativistic effective field theory, HB chiral perturbation theory rebuilds a power counting 
rule through simultaneous expansions in terms of 1 lm^ and external momentums. The pion-nucleon scattering has 
been investigated up to <9(/? 3 ) JslQll] and 0(p A ) | IT] with HB approach. Though the description of n-N scattering 



phase shift data is well described near the threshold region, the non-relativistic expansion encounters the problem of 
convergence in many cases flflll iljll. e. g. the scalar form factor of the nucleon does not converge in the region close 
to the two-pion threshold t - 4M% 



On the other side, in the framework of relativistic chiral theory, one may conclude that all the power-violating 
terms are polynomials and can thus be absorbed in the low energy constants from the effective lagrangian JT^Q- 
Hence the IR prescription and EOMS scheme are proposed to retain both correct power counting and covariance. 
Nevertheless, they are different in practice when removing chiral polynomials, the former subtracts all the so-called 
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infrared regular part of the loop integrals, which is always an infinite chiral polynomial of different order, while the 
latter only cancels the finite PCB terms. In Refs. |U 18] the pion-nucleon scattering amplitude is analyzed within 
IR prescription. In Ref. II 1611 . the 0(p 4 ) calculation was carried out and the analytic property of the amplitude was 
discussed. In Refs. 1171 Il8ll the 0(p 3 ) calculation result was used to fit the phase shift data. However, it has been 
shown that the pion-nucleon amplitude in IR prescription is scale-dependent fl6ll and suffers from an unphysical 
cut at u=0 fll8fl . Additionally, a huge violation of Goldberger-Treiman(GT) relation shows up [18], which queries 
the applicability of covariant baryon chiral perturbation theory. Hence these problems lead to the application of the 
EOMS scheme. 

The EOMS scheme provides a good solution to the PCB problem in the sense that it faithfully respects the analytic 
structure of the original amplitudes, e.g., see Ref. II 1911 . and being scale independent for merely making an additional 
subtraction of a polynomial of PCB terms with respect to the traditional substraction. A first attempt of the EOMS 
scheme on pion-nucleon scattering was made up to 0(p 3 ) level by Alarcon et al. i20ll . Remarkably, achievements 
in the description of violation of the GT relation 112 ill and the pion-nucleon sigma term Il22ll are obtained. However, 
as pointed in Ref. 1E3I1 . the size of the 0(p 3 ) contribution can be very large and comparable to those given by the 
lower-order terms even at very low energies above threshold. Thus, the applicability of EOMS scheme to describe 
the partial wave phase shift at <9(p 3 ) seems to be questionable. The authors of Refs. |2(| 23] solve the problem by 
explicitly including the contribution of A(1232) resonance. 

In this paper we extend the analysis up to 0(p 4 ), and settle down the convergence problem occurred in the 0(p i )- 
level analysis, even when the A(1232) resonance is absent. Compared with the 0(p 4 )-lR results in Ref. fl6ll . the 
pion-nucleon scattering amplitude presented here is the first analytic and complete 0(p ) result. Especially, we pay 
great attention to the substraction of PCB terms, such that the "threshold divergence" problem firstly pointed out 
by 12411 within the IR prescription never occurs. We perform fits to partial wave phase shift data and determine all 
the LECs involved. Besides, the leading order contribution from an explicit A(1232) resonance and partially-included 
A(1232) loop contribution are included to describe phase shift data up to energies just below the resonance region. 

A phenomenological discussion is also made based on the pion-nucleon scattering amplitude we obtain. We have 
mainly studied the Goldberger-Treiman (GT) relation and the pion-nucleon sigma term cr nN . The GT relation violation 
is a basic quantity to test the applicability of EOMS-B^PT to the pion-nucleon system. Hence we calculate it up to 
0{p y ). The prediction of the violation is in good agreement with other determinations II25[ 12611 . and its chiral series 



converges rapidly. The analysis on cr nN is important to understand the origin of the mass of the ordinary matter and can 
be useful for the study of the supersymmetric dark matter 1127[|28ll . Taking both pion-nucleon phase shift and the lattice 
QCD data into consideration, we give a reasonable prediction for the pion-nucleon sigma term: cr„ N -52 + 1 MeV 
from the fit without A(1232), and ct„n - 45 + 6 MeV from the fit with an explicit A(1232). The first one is smaller than 
the 0(p 2 ) result given by Ref. i20ll but larger than the recent 0(p 4 ) result in Ref. lE9ll . while the latter is in reasonable 
agreement with previous results found in the literature 11291 - 13 ill . 



2. Theoretical discussions on nN -* nN in EOMS scheme 



2.1. Kinematics and effective lagrangian 

In the isospin limit, the standard decomposition of the elastic n-N amplitude reads fil flill . 



Kn =xy\s*a'r + -[Tv.rjr- \ X N, r* = u{ P ', s') 



u(p, s) . 



(1) 



where p,q(p',q') denote the momenta of the incoming (outgoing) nucleons and pions, respectively, and a (a') stands for 
the isospin index of the incoming (outgoing) pion, see FigQ] T a >,T a are Pauli matrices and XN'iXn are the isospinors 



of the nucleons. For on-shell elastic scatterings, p 2 = p' 



, q 2 = q' 2 = M 2 and the Mandelstam variables s,t 



and u fulfill s + t + u 
D - A + vB with v = 



2m 2 N + 2M 2 „. Eq. ® can be written in an alternative form through the replacement of A by 



4m N ' 



T* = u(p', s') 



ir- 



2m ff 
2 



u(p, s) . 



(2) 




Figure 1: Kinematics of elastic 7T-N scattering. p,q (p',c/) denote the momenta of the incoming (outgoing) nucleons and pions, respectively, and a 
(a') stands for the isospin index of the incoming (outgoing) pion. 



Since the leading order contribution of A and B may cancel each other, one should better use D and B to perform 
the low energy expansion of the scattering amplitude when extracting the PCB terms. Our calculation of the 7T-N 
scattering amplitude up to 0{p A ) level demands the corresponding calculation of D ± and B ± up to 0{p 4 ) and 0(p 2 ), 
respectively. 

In chiral perturbation theory, each graph is assigned an overall chiral order D, which means the graph is of size 
(p/A) D , where one has the soft scale p <k A and A stands for a "high energy scale", i.e., the breakdown scale of the 
theory. For processes containing one baryon in the initial and final states, the chiral order for a given graph with L 
loops, V„ n-th order vertices, Nm meson propagators and Nb baryon propagators, is given by 

D = AL + J] nVn - 2N M - N B . (3) 

n 

The effective lagrangian relevant to the one-nucleon sector consists of n-N and purely mesonic lagrangian, 

r ff = r (l) + r (2) + r (3) + r (4) + ■ ■ ■ + r (2 > + r (4 > + ■■■ u) 

where the superscripts denote the chiral order. The lowest n-N lagrangian takes the standard form 

L^ = NUll)-m + ^g$-AN . (5) 

The nucleons are described by an iso-spin doublet as N - (n,p) T , and the covariant derivative acting on it is 
defined as = d M + F A ,, with 

= - [m t (5^ - ir^)u + u(_d^ - j'Z jU )m + ] , (6) 

where and r M are constructed from the external vector and axial vector currents as l M = v M - and r A , = v M + a^. 
The Goldstone bosons are collected in a 2 x 2 matrix-valued field u in the so-called exponential parametrization 

Uft = i [u' (d^ - ir^u - u(d^ - , u = ex P^"^y"j > ( 7 ) 

with t being the Pauli matrices. The parameters appearing in this lowest-order n-N lagrangian, m, F, and g are the 
bare values of the nucleon mass, the pion decay constant and the axial charge, respectively. 

For the complete form of -£S, •£ l and £^l, we refer to Ref. Il32ll . Here we only write down the terms which are 



3 



relevant to our calculation: 



4w = d<X+)NN - ■£^(u»u v )(ND ll D y N + h.c.) + ^<w%)AW - ^AyV Uy]N , (8) 



+ 7^3^. ID V , u A ]](D»D v D A + sym.) + h.c.) + i^-(\X-, u^W + h.c.) 
\2m i 2m 



+i — (o~ t ' v ([D / i, Uu\uy — Uu [D v , u A ])D' 1 + h.c.) 



.di 4 - d\s t 

— (IT \IL>A, i»iij«y — «uL 

8m 

+^Vr 5 Cr+>«, + ^yVcxW-lW , (9) 

4^ = N ( ei4 < - £^ {(hA,h\W v + h.c.) + ^ ({h^K^e + h.c.) 

+ ^[/^XK v - jrl ([h^Kp^D^ + h.c.) + e , 9 Cr + X« ■ u) 
~S (<*+>< l V M >'> Z> "' + h.c.) + ^Cr + >K, Mv ]^' v + e 22 [Z» p> [/y, <Jt+>]] 
"S (<^-V>^ V + h - c -) + ie i6 (u,[D^X-]) ~ [£>v,J-]]^ v + e38 <*+>Cr+> 

+ ^ Cr; -* 2 > - ^ (tr 2 > - Ct-> 2 + ixl) - Cr + > 2 )}^v , (10) 



■<xl -x 1 ) - e -f {(x 2 ) - ix-f + Or 2 ) - Cr + > 2 )}^v 

where the c,, <i, and e* are the low energy constants. The new symbols appearing here are defined as follows, 



X ± = u-\xu* + ux'u, 

h^y = [Dp, Uy] + [Dy, Up] , 

x± = x± - ^(x±) , 



= [D a D p ...D + sym) . (11) 



Hereof = M = diag(M 2 , M 2 ) and M is the bare pion mass. In the pure meson sector, the relevant terms of £^ and 
J$2 ar e given by 

4£ = ^4<«%><*+> + ^(/ 3 + / 4 )(y + > 2 , (12) 

and /3,Lt are low energy constants that will appear in our calculation, too. It is noticed that throughout this paper 
m, M, g, F represent the bare quantities for nucleon mass, pion mass, axial coupling constant and pion decay constant, 
respectively whereas m^, M„, gA, F„ the corresponding physical quantities. For the kinematic region close to the ttN 
threshold, one has 

a- . t M„ 

— «: 1 , — -« 1 , — « 1 , 
A 2 A 2 A 

or equivalently 

v . v B M„ 

— «; 1 , — <sc 1 , — •« 1 , 
A A A 

where cr — s - m 2 N , v = j^j*-, vg = ' ™J ■ Here the high energy scale A = {AnF„, m^, m&, m& - m^) with m& the mass 
of A(1232). Hence cr, t, M n (or v, v#, M„) are adopted as expansion parameters, and 

cr ~ O(p) , t ~ 0(p 2 ) , M„ ~ O(p) , v ~ O(p) , v B ~ 0(p 2 ) , m w ~ 0(p°) , m A ~ 0(p°) , m A - m N ~ 0(p°) . (13) 
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Figure 2: Tree graphs for 7T-N scattering up to 0(p ). The solid lines correspond to nucleons, while the dashed lines represent pions. The vertices 
with circled 2, 3 and 4 stem from XS, and -CLL respectively. The nucleon propagates with «!4 = m ~ 4c\M 2 - 2(8e38 + ens + eii6). Crossed 
diagrams for (a), (d) and (e) are not shown. 



2.2. Tree amplitudes 

We show all the tree graphs which contribute up to <9(p 4 ) in Fig. [2] according to the power counting rule given 
by Eq. (O, and list their contributions to D ± and B ± in | Appendix B| separately. The nucleon propagates with a mass 
parameter OT4 = m—Ac\M 2 -2(8e38 + en5 +en6)M 4 instead of m, so that the graphs with mass insertions, generated by 

Cl (x + ) in and e 3& (x + )(x + ) + e -f<xl ~X 2 -) ~ T (Ct 2 > - <X-> 2 + <xl> ~ <X + ) 2 ) in -C% in the nucleon propagators 
are automatically considered. For convenience, one can classify the tree graphs into two categories: Born-terms and 
contact terms. In Fig. [2] contributions from Born-term graphs, (a), (d) and (e), and their crossed diagrams can be 
summed and rewritten concisely in terms of the A and B functions as 



8 2 s 



- m 



2 



A ± = A(s) ± A(u), A{s) = — '-j (m 4 + m N ), 

gi s + 2m N ni4 + m 2 

B = B(s)±B(u), B(s) = -^ Z (14) 

4/* z s - m A A 

witlj|m4 = m - 4ciM 2 - 2(8e3g + ens + en6)M 4 and^2 = g + 2M 2 {2d\(, - c/ig). Meanwhile, the rest are contact term 
graphs without crossed diagrams, and the sum of them is 



16V 4 



4ciM 2 c 2 (l6m 2 N v 2 -t 2 ) c 3 (2M 2 -r) 4 , 2 

D+ = — W + g W + F 2 + J2 {<" K ~ *) + 2 ^ K -*)* + ^ 

D = ^ + J2 { 2 ^i +d2 + 2d ^ M * ~ + d ^ + 2d ^A ' 

+ 4(d u - d\s)vm N 1 2c 4 m N 8m N , , 2 \ 2] 

= J2 > = 2f2 H pi + ~pT r 17 \ 2 ^ ) + 2e ^ v } . (15) 

Note that in the graphs of 0(p 3 ) and 0(p 4 ) shown in Fig. [2] the bare constants can be replaced by the physical ones, 
since the distinction is beyond the accuracy of our calculation. Such replacements have been done in Eqs. [151 and 



2 As for the 0(p 4 ) effective lagrangian, we adopt the conventions of Ref. (32j], so ma, differs from the one in fl^l . where e^ L = -2(8^38 + ens + 
eil6>- 



Appendix B where the 0(p 3 ) and 0(p 4 ) contributions are expressed only by physical parameters. According to the 



discussion in Ref. Ill 111 , the terms proportional to (k = 19,20,21,22,35,36,37,38) only amount to quark mass 
corrections of c, (; = 1 , 2, 3,4), hence here we have already adopted in Eq. (fl~5T > the following combinations of LECs, 



ci - c\- 2M 2 (e 22 - 4e 38 ) , c 2 = c 2 + 8M 2 (e 20 + £35) , 

C3 = c 3 +4M 2 (2ei 9 -e 2 2 - <?36) » c 4 = c 4 + 4M 2 (2e 2 i - e 37 ) 



(16) 



2.5. Loop amplitudes 
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Figure 3: Loop graphs contributing to tt-N scattering at 0(p ). The nucleon propagates with 014 = m - \c\M - 2(8^38 + ens + fiii6). Crossed 
graphs for (a)-(i) and (n)-(s) are not shown. 



To carry out the calculation on jtN — > nN process up to 0(p 4 ) level, one ought to include 0(p 3 ) and 0{p A ) 
loop corrections with the corresponding Feynman diagrams shown in Fig. [3] and Fig |4] respectively. The 0(p 4 ) 
loop graphs in Fig [4] are simply obtained from the graphs in Fig. [3] involving even number of pions by replacing 
one of the O(p) vertices with corresponding 0(p 2 ) vertices. The results - in which no substractions have yet been 
performed - are displayed in Appendix C.2 and Appendix C.3 There the 0(p 3 ) loop results are listed explicitly for 
the sake of completeness. The definition of the loop functions is very similar to that in Ref. II 1611 . which is presented 
We have checked that our 0(p 3 ) loop results agree with those in Ref. 112 311 except a few terms due 



in 



Appendix C.l 



to the reason that we have chosen exponential parametrization instead of sigma parametrization for the pion fieldsrTi 



J All the terms proportional to f~g in Ref. fl^l should be reversed in sign, as it is firstly pointed out by Ref. flUl . 

4 The physical pion-nucleon scattering amplitudes up to 0(/? 3 ) here and in Ref. I23I1 are the same and are independent of parametrization. 
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Figure 4: Loop graphs contributing to ^-N scattering at 0(p 4 ). The vertex with circled '2' stems from The nucleon propagates with 

nt4 = m - 4ciM 2 - 2(8^38 + ens + eii6)- Crossed graphs are not shown. 



To our knowledge, the Q(p 4 ) loop contributions shown here are the first analytic and complete calculation results, 
and one can consult in Ref. Il 111 for the 0(p 4 ) results within heavy baryon^-PT and in Ref. [ 16] for the 0(p 4 ) results 



within infrared regularization. In Ref. I116ll special technical simplification of calculation is adopted, here we refer the 
readers to section 6 there for details. 

It must be emphasized that the parameter m in both 0(p 3 ) and 0(p 4 ) loop results should be understood as m\ = 



m - Ac\M 2 - 2(8e38 + C115 + e\i^)M 4 for the following reasons f 33] : 

• All loop diagrams with contact interaction insertions in the nucleon propagators are summed up automatically; 

• When performing renormalization, one can directly set m.4 = in the one-loop results up to 0(p 4 ) level, since 
corrections are at least of two-loop order (the lowest chiral order of two-loop contributions is naively 0(p 5 )). 

2.4. EOMS scheme and PCB terms 

Since the nucleon mass my is nonzero in the chiral limit, the necessary power counting rule for an effective theory 
breaks down, namely PCB problem occurs. To remedy the PCB problem we adopt the EOMS scheme proposed by 
T. Fuchs et al. JH, which suggests performing renormalization in two steps: the firstly traditional MS-1 substraction 
to cancel the ultraviolet divergencies and then EOMS substraction to remove the PCB terms. The EOMS substraction 
is remarkable in the sense that the renormalized nN — > nN amplitude will possess good analytic and correct power 
counting properties since the PCB terms are polynomials of quark masses and momenta and are absorbed in the LECs 
eventually. Especially, as proved by Bee her and Leutwyler |]6j], the PCB terms stem from the regular part of the loop 
integrals, which allows us a simple way to obtain the PCB terms if we have known all the regular parts of the loop 



integrals needed - these are shown in Appendix D 

Taking the 0(p 3 ) loop amplitude for example, one firstly changes the amplitude in A, B form to D, B function^ 
and reduces them to expressions only containing scalar one-loop integrals. Then those scalar one-loop integrals are 
substituted by their regular parts to a given order, and a chiral expansion in terms of small quantities like M,t,cr = 
s - mr is performecj^. Finally, for the total 0(p 3 ) loop amplitude, the series whose chiral order are lower than 0(p 3 ) 
are regarded as PCB terms which read 



5 In Ref. the fact that the leading contribution from A and B cancels, while not for D and B, is pointed out. 

6 One can also chose cr = s - s,i,, with s,i, = (m + M) 2 , as expanding parameter like Ref. |24|| . here we follow Ref. Q]. 



£>C3H = J__ { 6 ^2 m 2 M2(r 2 + 2(J .4 + g 4 J 2m 4 ( 1Mf » _ ?M 2 f + ^2 (_ ?M 2 + 3^ ^2 + ^4] j > 



64F 4 

D ^ = K ( _2m2 +t+2 °)- 2mZ ( 2m2 - K - f + 2 -)} • 



BpcB = K -* + *r). = ^f^2 K + g2 K (" 5M2 + f ) + • (17) 

The same procedure can be taken to extract the PCB terms of the total 0(p ) loop amplitude, 

D% = n52 l^ 2(r 3 (864cig 2 m 2 M 2 cr 3 + [l6c 4 - (9c 2 - 216c 3 + 272c 4 )g 2 ] m\t - 2M 2 )<? 

-(9c 2 + 32c 3 + 32c 4 )g 2 m 4 (f - 2M 2 ) 2 (2M 2 - t + (f) - 2(9c 2 + 32c 3 + 32c 4 ) (? 2 mV 4 
+4 [2c 4 - (9c 2 + 16c 3 + 14c 4 )g 2 ] (f - 2M 2 + cr)cr 4 } , 

O^s = 2304^V 3 ( 32(2C2 + 1?C3 ~ l9c ^8 2 m 2 M 2 (2M 2 - f)tr 2 

-4 [l44ci - 2c 3 - c 4 + 4(18ci + c 2 - c 3 - 7c 4 )g 2 ] fcr 4 

+8 [(72 Cl - 2c 2 - 6c 3 - 3c 4 ) + (36cj - 2c 2 + 9c 3 - 33c 4 )g 2 ] (f - 2M 2 )cr 4 

-(9c 2 + 32c 3 + 32c 4 )g 2 m 2 (f - 2M 2 ) [2m 2 (2M 2 - f)(2M 2 - t + cr) + cr 2 (4m 2 + cr)] 

+ ^2 [(34c 2 + 30c 3 - 3c 4 ) - (4c 2 + 15c 4 )^ 2 ] <x 6 1 , 



B Prn = t~i — T (24c 4 cr 4 + g 2 [32(2c 2 + 17c 3 - 19c 4 );7? 2 M 2 cr 2 + (67c, - 56c 3 + 96c 4 )cr 4 l 

576/ 4 7r z cr- i 1 L J 

+2(9c 2 + 32c 3 + 32c 4 )g 2 m 4 [4M 4 + t 2 - tcr + 2cr 2 + M 2 (-4f + 2cr)]} , 

3 

fi PCB = 576 J 4?rV3 {[ 9c 2 + 32c 3 + !6c 4 - 2(9c 2 + 16c 3 - 28c 4 )g 2 ] a 3 

+2g 2 m 2 (9c 2 + 32c 3 + 32c 4 )(2M 2 - f)(2M 2 - / + oj\ . (18) 

Those PCB terms will be substracted, namely absorbed by redefinition of the LECs, when performing the EOMS 
renormalization of the nN — > ttN amplitude in section [231 

Before ending this subsection it may be worthwhile to mention that the PCB terms should be prevented from 
divergences induced by prefactors of the type 

1 1 

and — =— =- , (19) 



A{s,m 2 ,M 2 ) /i(s,m 2 ,M 2 ) + st' 



respectively, and so should the EOMS-renormalized amplitude be. The obstacle is first noted by Ref. [24], that 
the numerial analysis of the IR-renormalized amplitude encounters divergences at threshold s t h = (m + M) 2 and at 
t — -A(s,m 2 , M 2 )/s. Nevertheless, both the EOMS- and IR-renormalized amplitudes should possess good analytic 
properties at those s and t values, namely singularities caused by ([T9T l are canceled by the numerators of the amplitudes. 
It can be easily seen from HB(s) { '\i = 1, ■ ■ ■ ,6) and H { ^(j =1,2) in | Appendix C.l that the prefactors are actually 
introduced by the standard Passarino-Veltman decomposition \34l of tensor integrals, which can be avoided by the 
new approach developed in Ref. l35ll . 

Taking the tensor integral for example, in Passarino-Veltman approach Hg is decomposed into 

where the expressions for and refer to Eq. (IC.5b and Eq. (1C.61 . On the other hand, following the approach in 
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Ref. 13511 . H^ B is now decomposed into the new form 



H 2 



n r a"- 

7 J (2^r 

n r d d - 

7 J (2* 



k 1 



)d + 2 [ M 2 _ £2] [ m 2 _ (p _ fc)2]2 [ OT _ (S - fe)2] 
d+2 k 1 



)rf+2 [ M 2 _ £2] [ m 2 _ (p _ fc)2] [ OT _ (£ _ fc)2]2 



where //[ and H 2 are already scalar loop integrals in dimension 6 momentum space, in other words, coefficient like 
1 /A(s, m 2 , M 2 ) in standard Passarino- Veltman decomposition never occurs. Hence the threshold divergence introduced 
by 1/A(s, m 2 , M 2 ) disappears, the same conclusion holds for the divergence at t — -A(s, m 2 , M 2 )/s when considering 



The new approach enables us to reduce the tensor integrals in the amplitude to scalar integrals defined in higher 
dimension momentum space, without confusion such as the divergence at threshold. If further regular parts of the 



scalar integrals are known, the PCB terms can be obtained in a new way. In Appendix D the method proposed by 
Refs. fl pOll is adopted to calculate the regular parts of the scalar integrals in dimension 4 space, those of the scalar 
integrals in higher dimension space can also be calculated term by term using the same method. With the aid of the 
regular parts, the PCB terms for Hi and H 2 can be easily obtained, which are 3 , 1 2 an d W^n? ' respectively- Using 
the relations = \{H X + H 2 ) and Hf } = \{H X - H 2 ), one then gets the PCB terms for H^ and Hf } in this new 
way, which are 30 \_ m i and respectively. On the other hand, the PCB terms for H^ and H ( ^ can also be obtained in 
the usual way adopted in the paper. One firstly replaces the scalar integrals in Eqs. IC.5l and lC.6l bv their regular parts 
shown in | Appendix D| then expands and H^' in terms of cr, t, M„. From the expanded expressions of and 
H^ 2 \ one finds that the PCB terms for H^ and H {2 ' are 31 \_ m i an d 0, which are the same as the expressions obtained 



through the new way above. In this sense, the approach developed in Ref. 135I1 provides us a new way to obtain the 
PCB terms. 

2.5. Renormalization 

As an example, we will firstly show the renormalization of nucleon mass as well as axial-vector coupling g& 
to interpret the essence of EOMS scheme. Noticeably the expressions of and g& are also needed for replacing the 
correspondin g ba re quantities in the tree amplitude when performing numerical fits. Part of the results are already 



given in Ref. IB6I1 



2.5.1. Nucleon mass and wave-function renormalization constant 



(a) (b) (c) 

Figure 5: One-loop Feynman diagrams contributing to the self-energy of the nucleon up to 0(p 4 ). The vertex with circled '2' stems from X!S. 

As a simple example to illustrate the EOMS method, we evaluate the nucleon physical mass up to 0(p 4 ). The 
0(p 4 ) result of m^ in the EOMS scheme can also be found in Refs. H, 37]. The one-loop Feynman diagrams are 



depicted in Fig. [5] A straightforward calculation leads to the primitive expression for the nucleon mass, 

3'V U „2„, 20 . 3M 2 



m N 



m - 4ci M 2 - 2e m M 4 - {A N - M 2 H(m 2 )} + — {(2ci - c 3 ) - ^} A* 



2„2 



3c 1 M I g 
F 2 



[-2 [A N - M 2 H(m 2 )] + 4m 2 [^(0) - M 2 H A (0)]) , 
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where e m = 8e38 + ens + eii6- The second and third terms are tree contributions stemming from 0{p 2 ) term involving 
Ci and 0(p 4 ) term related to e^, em and ene- The term with the first, second and third brace bracket represents 
the loop contribution from (a), (b) and (c) in Fig. [5] respectively. Definitions of all loop functions appeared here 
follow Appendix C.l One can perform different renormalization schemes on the above expression of m^, e.g. IR 
prescription, etc.. However, we proceed with the EOMS remormalization by firstly carrying out traditional MS - 1 
substraction, which gives 



m N 



= m r - Ac\M 2 - 2e' m M 4 - |?£ [A N - M 2 H(m 2 )} + ^ {(2c, - c 3 ) - °j] A„ 



2„2 



3ciM z g 
F 2 



[-2 [A N - M 2 H(m 2 )] + Am 2 [J N (0) - M 2 H A (0)]} , (20) 



The bar over the loop function denotes the finite part of it, and the LEC with a subscript r means that it is a MS - 1 
quantitjQ. The MS - 1 substraction does nothing but shifts the divergencies in loop functions to the bare mass and 
LECs: 

3m?g 2 R 
m — m — , 

32F 2 

Cl = Cl+ 128^ (1_12Cim) ' 
3R 

em = e m + mF 2 n 2 h 8 ^ 1 + ^ + ^2 + 4c 3 ] , 

where R = -r? + y^ — 1 — In An. Since m# is scale-independent , we now take the renormalization scale p — m in 
Eq. ( 1C.4I ) for simplicity, and therefore A# = 0. If the loop functions are replaced by their regular parts, one naively 
finds that the term im ^. 2 8 H(m 2 ) in Eq. d20b should be 0(p 3 ), but actually contributes an 0(p 2 ) PCB term ^^tjt, and 

the same thing happens for the last term where a 0(p 2 ) PCB term 3c ' 4 ^. 2 ^ 2 8 occurs. Since they are polynomials, they 
can be absorbed by the LEC c\, 



3mg 2 

C ~ 1=C ^128^ (1+8Cim) ' 
Finally, we get the expression for in EOMS scheme, which takes the following form 



3mM 2 e 2 — , 3M 2 ( C7 , 

m - Ac x M 2 - 2e m M 4 + ^ H(m 2 ) + — {(2c, - c 3 ) - -j \ A n 



K « (m!) + W [« 0) - „^ (0)] ) - fgv _ 3^ m 

e r m , namely they are unaffected by the PCB terms. 
We note that one can also carry out the mass renormalization by replacing m by = m — 4 Cl M 2 - 2e,„M 4 m 
the nucleon propagator. In this case, the graph (c) in Fig. |5]is absent and automatically included in graph (a), and the 
result is 

m N = m - 4c,M 2 - 2e,„M A - ^f- {A^ - M 2 H(m 2 A )) + {(2d - c 3 ) - °j } A. , (22) 
while the wave-function renormalization constant of nucleon reads 

Z N = 1 - % L - AmlM^His)) -p-^. (23) 
AF 2 { ds J f=mN F 2 m 4 d 

Hereafter, the 1114 related to loop contributions is always taken as m for short. Instead of Eq. (fJTJ, Eqs. (l22l i and (l2Jt are 
adopted for the renormalization of the n-N scattering amplitude. Throughout this paper, we use this way to simplify 
our calculation for the reasons discussed in section [ 
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(a) 



(b) 



[ 1 (c 



(d) 



(e) 



(f) 



(g) 



Figure 6: Diagrams contributing to the nucleon axial form factor Ga(0 up to 0(p ). The wavy line with a circled cross at the end stands for the 
axial-vector current. The vertex with circled '2' and '3' stems from xS and ^52, respectively. 



2.5.2. Axial-vector coupling constant 

The axial-vector current A^ "(0) between one-nucleon states can be written as 



(N(p')\A»- a (0)\N(p)) = u(p') 



2 "> 1 

GU(<7 )Y I 75 + G P {q-)—y 5 
2m 



t" 



-u(p), 



(24) 



where = p'^-Pf, and a is isospin index. GAiq 2 ) is called the axial form factor and Gp(q 2 ) is the induced pseudoscalar 
form factor. The axial-vector coupling constant gA is defined as 



g A = G A (q 2 = 0). 

Up to the 0(p 4 ) diagrams in Fig.|6]are needed. A straightforward calculation leads to gA'- 



(25) 



g 3 m 2 g(4-g 2 ) A g(2 + g 2 ) A g\2m 2 + M 2 ) 

■N ' 



g A = g + 4d i6 M 2 - o ° ^ - + y A;v - °\^° ' K + - v "~ A ^ — ' Jn(0) 



32F 2 tt 2 2F 2 
g(8-g 2 )M 2 rT/ . # 2 M 4 



2F 2 
3g 3 ra 2 M 2 5 



4jF 2 



4F 2 "("V^W ^ & 



8c 4 mg 



(A, - M 2 H(m 2 )) 



2g 
mF 2 



c 2 U— A„ + —Aw - M 2 i/ (2) (m 2 )) - 4(c 3 + c A )m 2 H (2 \m 2 ) 
\ d d j 



(26) 



which agrees with Refs. J38, 39 1 Here, the EOMS renormalization procedure is similar to that for m^. The bare g 
will be redefined as 

(9c 2 + 32c 3 + 32c 4 ) ,f = g+ t|^(2 - g 2 ) - J^(3c 2 + 8c 3 - 40c 4 ) , (27) 



6 6 16F 2 tt 2 576F 2 tt 2 
whereas the redefinition of the LEC c/ie reads 



16F 2 7T 2 



(c 2 + 18c 3 - 18c 4 ) , d[ 6 =d l6 + 2 [3g(l - g 2 ) - mg(c 2 + 6c 3 - 18c 4 )] . 



96F 2 7r 2 



H6 - a i 6 



288Ftt 2 



\92F 2 n 



(28) 



The final expression for gA is lengthy but rather straightforward to get with the help of Eqs. (l27b and d28l ), so we do 
not present it here explicitly. 



'Hereafter we denote the MS — 1 and EOMS renormalized LECs with a superscript r(eg. Cj ) and overhead tilde (eg. c~i) respectively. 
8 Though in Ref. (3^1 the mass insertion graphs are calculated directly, the comparison between our result and the one there is easy. 
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2.5.3. Full nN — > nN amplitude 

In order to present a full nN — * /rW amplitude, we also need the formulae for M n and F„, and the corresponding Z 
factor for the pion . To 0{p 4 ) level, these read 



Mi 



M 



1 + 2h^r + 

f 2 2/2 



F„ = F 



M 2 



1 



1 +4— r - -rA 



1 + 



/ 2 



f 2 



/ 2 



(29) 



All of them do not contain PCB terms from loop integrals, and hence can be treated traditionally. 

Since all the necessary preparations are completed, we proceed with the renormalization of nN — > nN amplitude. 
Unlike the renormalization of and gA, it is hard to visualize the procedure of nN — » nN amplitude renormalization 
for its extremely lengthy expression. However, the essence is the same, that is to carry out renormalization procedure 
in two steps: MS - 1 renormalization and EOMS renormalization. Corresponding to the MS - 1 renormalization, 
those LECs appeared in the tree amplitudes are demanded to cancel the ultraviolet divergences and yield the so-called 
MS - 1 renormalized LECs, 



c'M) 



■ypnR 
l6n 2 F 2 ' 



l6n 2 F 2 ' 



el(M) = e k 



\6n 2 F 2 m 



where details of y^,y^,yl can be found in Appendix E To absorb the PCB terms (fTTt and (fT8T l, c[(//) and d r j(p) are 
further redefined as 



C; = C: - 



dj = d r j- 



6 c .m 

l6n 2 F 2 ' " J ~ " } \6n 2 F 2 

e r Ap), since the chiral order of PCB terms are lower than 0(p 4 ). 



whereas e r k (p) remain the same, ie. e\ = e r k (p), since the chiral order of PCB terms are lower than 0(p ). Also 
61 and dj are specified in Appendix E So far, we have already completed the renormalization of the nN — > nN 
amplitude in the EOMS scheme, the main feature of this method is characterized by additional EOMS substractions, 
which distinguishes EOMS scheme from other prescriptions like IR and HB. We observe that an amplitude in EOMS 
scheme differs from the full covariant amplitude only by a polynomial of small quantities and hence owns the same 
analytical structure but possesses correct power counting . The validity of the nN — » nN description in EOMS scheme 
will also be judged by numerical fits to existing experimental data. 

2.6. Partial wave expansions 

We choose to perform fits to the partial wave phase shift data. The iso-spin decomposed amplitudes for nN 
scattering are 

T /= 5 = T + +2T~, 

T !=3 2 = T + -T~. (30) 

The final partial wave amplitudes with isospin / , orbital momentum 6, and total angular momentum / = t± \ (denoted 
by l± concisely) take the form (lj], 

lo7rys 1 L v ' J 



where 



+(E P - m N ) \-A' (±l (s) + ( aAs + m N ) fi' ±1 (s)]} , 

A' e (s) = j A I (s,t)P e (cos6)dco&e , 

B' e (s) = J B'(s,t)P e (cos6)dcos6 . 
12 



(31) 



Here E p = J+ ^ l ^- M " and 6 are the nucleon energy and scattering angle in center-of-mass system, respectively. fV(cos 0) 
are the conventional Legendre polynomials. The angular variable cos 9 relates to the Mandelstam variables via cos = 
1 + — H ' M iy with A(a, b, c) = a 2 +b 2 +c 2 -2ab-2bc-2acbeingthe Kallen function. As a straightforward consequence 

of unitarity of S matrix, one can further express the partial wave amplitudes by the phase shift S 1 . , 

fL(s) = [exp(2wSj ± (j)) - l] , (32) 

where p is the 3-momentum of nucleon in the center-of-mass frame. Since the phase shift is real for elastic scattering, 
we follow Ref. ifioll to related it with our perturbative computation of f' ± (s) via 

4(*)=arctan{|^Re// ± 0s)} . (33) 



3. Phenomenological and numerical studies 

In this section, we firstly perform fits to partial wave phase shift data near threshold to pin down the free LECs. In 
order to describe the partial wave phase shifts up to a higher energy region, we include explicitly the leading A(1232) 
Born-term contribution and partially-included A(1232) loop contribution. The contribution to the LECs from the 
A(1232) resonance is also considered. We proceed with discussing the convergence of the chiral expansion of the 
resulting partial wave phase shift. The improvement of the fourth-order calculation compared with the third-order is 
shown. Finally, the deviation (Act) of Goldberger-Treiman relation and the pion-nucleon cr term a^N are discussed. 
The 0{p 3 ) analyses are also included for the sake of comparison with the previous literature. 



3.1. Partial wave phase shift 

To begin with, we firstly fit the partial waves at the 0(p 3 ) level. We denote this fit by "Fit \-0(p 3 y\ As input we 
use the phase shift data from Ref. [40], namely the current solution of George Washington University (GWU) group. 



Since the GWU group does not give data errors, we assign them with the method of Ref. 11811 . 



r(<5) = ^e] + e 2 5 2 , (34) 



with the systematic error e s = 0.1° and the relative error e r = 2%. Throughout the numerical analyses, we employ 
g A = 1.267, F„ = 92.4 MeV, m N = 939 MeV, M n = 139 MeV, and the renormalization scale p = m N . There are 9 
free LECs (or combinations of LECs) in total: ci_4, d\ + di, di, d$, du - d\$, if is. All of them can be pinned down 
by fitting two S- and four P- partial waves. The fitting range is from threhold (1.078 GeV) up to 1.130 GeV in yfs, 
and the interval between two data points is 4MeV. The 2nd column in Table [T]collects our fit results at 0(p 3 ) level. In 
column 3 of TableQ] we have also listed the results from Refs. 112 ill for comparison. We see that, in general, our fit 
results at 0(p 3 ) level are in good agreement with those in Refs. 112 ill , except for the ds parameter. Especially, the d\%, 
related to Act, is nearly the same. 

The fourth order analysis of n-N scattering is denoted by "Fit I (a)-<9(/? 4 )" in Table 2. There are 14 free LECs, 
which are four dimension two LECs: c\, q, £3 ,£4 , five dimension three LECs: d\ + di, d^, ds, d\\ - dis, d\%, 
and five dimension four LECs: e^, e\s, eie, en, e^. Unlike the 0(p 3 ) fit, d\% is now fixed at its 0(p 3 ) fitted value, 
according to the discussion of Act below in Section [3~4l The 0(p 4 ) fit is performed up to 1.1 3GeV too, and the results 



are shown in column 2 of Table [2] Also, we have taken the results of HB^PT from Ref. Ill ill for comparison. Our 



results show improvements compared to Ref. ill ill . Firstly, from Table [2] one can observe that the c,, dj and e^ are 



mostly of natural size in EOMS scheme, but in HB results, especially Fit 2 in Table 1 of Ref. Ml 111, some of the 
come out fairly large. Secondly, our results seems to be more self-consistent. The c,- change a lot when extending the 
0(p 3 ) analysis to the 0(p 4 ) analysis in Ref. fllll . while our c, change much more acceptably. 

We plot both Q(p 3 ) and Q(p 4 ) fits in Fig. [7] Though fits are performed up t/s — 1.13 GeV, we plot up to 1.16 GeV. 
The conclusions made in Ref. Ill 111 still hold: the P33 wave is slightly improved compared to the 0(p 3 ) calculation, 



and the fn partial wave are somewhat off above 1.14 GeV. 

Note that the effect of unitarity is automatically included through the phase shift formula Eq.[33] which is discussed 



in Appendix F 



13 



LEC 


Fit l-Oip*) 


WI08 [21J 


Fit ll-Oip 3 ) 


WI08 [23J 


C\ 


-1.39 ±0.07 


-1.50 ±0.06 


-0.81 ±0.03 


-1.00 ±0.04 


Cl 


4.01 ± 0.09 


3.74 ± 0.09 


1.46 ±0.09 


1.01 ±0.04 


Ci 


-6.61 ±0.08 


-6.63 ± 0.08 


-3.10±0.12 


-3.04 ±0.02 


C\ 


3.92 ± 0.04 


3.68 ± 0.05 


2.35 ± 0.06 


2.02 ±0.01 


d\ + di 


4.40 ± 0.54 


3.67 ± 0.54 


0.79 ± 0.09 


0.15 ±0.20 


dj, 


-3.02 ±0.51 


-2.63 ±0.51 


-0.47 ± 0.05 


-0.23 ± 0.27 


ds 


-0.62 ±0.13 


-0.07 ±0.13 


-0.17 ±0.04 


0.47 ± 0.07 


di4 - dis 


-7.15 ± 1.06 


-6.80 ± 1.07 


-0.90 ±0.15 


-0.5 ± 0.5 


di& 


-0.56 ± 1.42 


-0.50 ± 1.43 


-0.91 ±0.25 


-0.2 ±0.8 


h A 




2.82 ± 0.04 


2.87 ± 0.04 


Xd.o.f 


0.20 


0.22 


0.35 


0.23 



Table 1: LECs given by fit up to <5(p 3 ). Fit I is performed up to 1.13GeV, while Fit II up to 1.20GeV including the explicit A(1232) contribution. 
For comparison, we provide the results from Il2lll23ll . The c, and di have units GeV~' and GeV~ 2 respectively, and h\ is dimensionless. In fit II, 
the results correspond to c' [ and d'j instead of c\ and dj, respectively. 



3.2. Contribution of A(1232) 

In this subsection, the effect of A(1232) is explicitly included to describe the partial wave shift up to 1.20 GeV. 
Pascalutsa et al. discussed how to treat the A(1232) as an explicit degree of freedom in covariant baryon chiral 



perturbation theory in Refs. H41N43H . The description of A(1232) is subtle, because the conventional Rarita-Schwinger 
representation is a field with 16 components while only 8 of them are physical. However we adopt the consistent 
formulation here l42ll . Additionanlly, we follow the so-called 5-counting rule l4lll which assigns an extra factional 
suppression of 0(p 1 ^ 2 ) to the propagator of A(1232). Up to 0(p 4 ) level, there are three typical A(1232)-included 
Feynman diagrams of different order: Born-term of 0(p 3 ^ 2 ) and 0(p 5 ^ 2 ), loop graphs of 0(p 7 ^ 2 ). Refs. |2^, 23 1 
remarked that the contribution of Born-term of 0(p 5/2 ) is negligible. It is rather complicated to evaluate loop graphs 
of 0(p 7 ^ 2 ) in the EOMS scheme, since the loop diagrams involving both propagator of nucleon and A( 1 232) will cause 
much more subtle PCB effects due to the heavy masses and m&. So throughout this paper we consider the leading 



Born-term contribution of A( 1 232), whose expression can be found in Appendix A. 1 together with partially-included 
A loop contribution illustrated in Appendix A. 2 The complete calculation with A(1232) up to 0(p 7 ^ 2 ) is left as an 



open question. It is important to mention that the effect of the A(1232) width is considered through the phase shift 
formula Eq. [33] which is amply discussed in Appendix G Likewise, we will have the operators from Eqs. (ISl-dTOt 



but with couplings different from those in the A-less effective field theory. We will mark the analogous coupling of 
the theory with A(1232) with a prime, e.g., c,- — > c' r 

Corresponding to the two different fits performed in Sec. 13.11 we perform another two fits, "Fit R-0(p 3 )" and 
"Fit II(a)-(9(p 4 )", which explicitly include the A(1232) contribution. The results are shown in the 4th column of 
Tableland Table|2j respectively. Taking into consideration the A(1232) contribution to the LECs, the d\% is set free 
in both fits here. The leading Born-term contribution of A(1232) is characterized by the NA coupling h A . The value 
of h A = 2.90 is determined from the Breit-Wigner width T A = 118MeV. In "Fit II(a)-0(/)", we fix h A = 2.90. 
However, for "Fit ll-0(p 3 )", h A is released as a free parameter, and its fitted value is 2.82. In Table [TJ our result is 



found mostly compatible with those of Ref. 12311 . We plot the A-included 0(p 3 ) and 0(p 4 ) fits together in Figj8]for 
the convenience of comparison. We find that, both 0(p 3 ) and 0(p 4 ) calculations with A(1232) contribution give a 
reasonable description to data and the 0(p 4 ) calculation improves the fit quality. 

Note that we denote the 0(p 4 ) fits with the notations: Fit I is A-less and Fit II includes A; Fit (a) and Fit (b) are 
performed only with pion-nucleon phase shift data whereas Fit (b) has the c, in c,- fixed, and Fit (c) is performed with 
both phase shift data and lattice QCD data for m^. 

3.3. Convergence properties of partial wave phase shifts 

We have made chiral expansions up to 0(p 4 ). It is necessary, at this stage, to check the convergence property of 
the chiral amplitudes. However, in the fits of Table [2]the 0(p 2 ) parameters c, mix with 0(p 4 ) parameters e,, so it is 
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LEC 


Fit l(a)-0(p 4 ) 


hb^pt yjj 


Fit II(a)-0(p 4 ) 


Fit l(c)-0(p 4 ) 


Fit ll(c)-0(p 4 ) 


C\ 


-1.08 ±0.06 


(-3.31,-0.27) 


-1.03 ±0.03 


-1.09 ±0.08 


-0.95 ± 0.05 


Cl 


2.78 ±0.11 


(0.13,3.29) 


0.50 ± 0.04 


2.44 ± 0.05 


0.10 ±0.06 


Ci 


-5.26 ±0.14 


(-10.37,-1.44) 


-3.17 ±0.05 


-5.05 ±0.22 


-2.64 ±0.08 


C4 


2.43 ±0.19 


(2.80,3.53) 


0.79 ± 0.03 


2.43 ±0.19 


0.80 ±0.03 


d\ + do 


6.29 ±0.12 


(4.45,5.68) 


2.99 ± 0.05 


6.18 ± 0.11 


2.93 ± 0.05 


di 


-6.87 ±0.16 


(-4.91,-2.96) 


-5.04 ±0.05 


-6.87 ±0.15 


-4.90 ± 0.04 


ds 


0.51 ±0.11 


(-0.95, -0.09) 


1.32 ±0.04 


0.55 ±0.11 


1.24 ±0.03 


du - dis 


-12.09 ±0.24 


(-11.14, -7.02) 


-5.61 ±0.09 


-11.94 ±0.23 


-5.58 ±0.09 


di& 


-0.56* 


(-1.53,-0.85) 


1.14 ±0.20 


-0.56* 


1.64 ±0.17 


eu 


3.69 ±0.36 


(-4.68,7.83) 


-4.53 ± 0.09 


-1.80 ±0.33 


-8.22 ±0.08 


<?15 


-14.99 ±0.55 


(-18.41,9.72) 


5.05 ±0.13 


-5.41 ±0.57 


10.52 ±0.12 


Cl6 


7.35 ±0.35 


(6.42,7.79) 


-0.31 ±0.07 


4.34 ±0.28 


-1.50 ±0.05 


en 


-2.29 ± 1.34 


(-17.79,14.88) 


16.98 ±0.15 


-2.23 ± 1.42 


15.70 ±0.15 


ei8 


6.07 ± 1.18 


(-9.15, 19.66) 


-10.99 ±0.12 


6.00 ± 1.26 


-9.87 ±0.12 


h A 






2.90* 




2.90* 


e\ 








15.48 ±0.30 


16.70 ±0.27 


m 








0.88 ± 0.02 


0.89 ± 0.03 


Xd.o.f 


0.04 


(0.008,0.44) 


0.23 


0.51 


0.36 



Table 2: LECs given by fit up to 0(p 4 ). Fit 1(a) and Fit 11(a) are performed with phase shift data, while Fit 1(c) and Fit 11(c) with both phase 
shift data and QCD lattice data (see section [331 . Fit 1(a) and Fit 1(c) are performed up to 1.13GeV, while Fit 11(a) and Fit 11(c) up to 1.20GeV 
including the explicit A(1232) contribution. In fit 11(a) and fit 11(c), the results correspond to c' r d'. and e' k instead of c,, dj and ek, respectively. For 
comparison, we provide the results from (Till . The e;, dj and have, respectively, units of GeV -1 , GeV~ 2 and GeV~ 3 , and Iia is dimensionless. 
The * denotes an input quantity. 



not suitable for testing the convergence of the chiral expansion. To overcome the problem, here we follow the strategy 
of Ref. Ultl to redo the fits. That is, we fix cj_4 in ci_4 of Eq. (TTST i with their corresponding fit values at 0(p 3 ) level 
given in Table [TJ In other words, we perform fits with four dimension-4 combinations: en - 4e^ (in ci), e^o + <?35 
(in ci), 2e\g - en — e^ (in C3 ), 2e2i - ^37 (in q), instead of ci_4. For clarity, we will denote the modified 0(p ) fits 
discussed here by "Fit I(b)-<9(p 4 )" and "Fit II(b)-(9(/? 4 )", respectively. In this case, the contributions from different 
orders are separated, so we can study the convergence of the amplitude. The resulting values for the LECs are shown 
inTable|3] 

Comparing "Fit l(b)-0(p 4 )" results with those from Ref. lfTlll . which are summed as intervals listed in the third 
column of Table [3] it is found, however, that most of our fitted dj and do not locate inside the intervals. The main 
reason might be that our primitive values for C2-4 as input are not in the corresponding intervals (see Table [3}, which 
cause the incomparability since a small variation of c, may lead to big changes of the higher order LECs, dj and e^, 
though both our fit and that of Ref. [11] maintain a good convergence property. 

The convergence can be visualized by plotting contributions from O(p), 0(p 2 ), 0(p 3 ), 0(p 4 ) separately, and the 
sum of them in Fig. [9] Note that we plot up to 1.20GeV, though fit only up to 1.13 GeV. One can observe that the 
0(p 4 ) contributions (cyan dashed-dotted lines in Fig. [9j» are in general small for all the partial waves below 1.13 GeV. 
The 0(p 3 ) contributions (magenta dotted lines) are mostly larger than O(p) contributions (green dashed lines), with 
the exception of the S n and P33 partial waves. However, there exist cancellations between the 0{p 2 ) (blue shorted 
dashed lines) and 0(p 3 ) contributions. The red solid lines represent the total contribution up to 0(p 4 ). They describe 
the existing partial wave data below 1.13 GeV very well. After all, the convergence property of the fourth order 
calculation is reasonable, while the third-order calculation is not satisfactory as pointed out by Ref. i23ll . 

In Fig. [10] we include the A contribution and plot the contribution from A(1232) together with the contributions 
from O(p), 0(p 2 ), 0(p 3 ), 0(p 4 ), and the sum of them. From the yellow short dash-dotted lines in Fig. [10] we can 
observe that A(1232) mainly contributes to f^-wave while the contribution for other channels is non-zero but very 
small. On the other hand, the chiral series are in general well convergent near threshold almost for all the partial 
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Figure 7: (Color online) Fit up to 1.13 GeV. The fourth- and third-order fits are presented by the blue solid and green dash lines respectively. 
Results of Fit l(a)-0(p 4 ) in tableland Fit l-0(p 3 ) in table[T]are adopted for plotting. 



waves. However when increasing the energy far above the threshold, the convergence becomes worse. Especially, we 
can see from the Fig. [10]that the higher chiral order contributions grow much more rapidly than the lower chiral order 
contributions as the energy increases. This indicates that the chiral perturbation expansion breaks down in the large 
energy region and stops being valid. 

3.4. Goldberger-Treiman relation 



The Goldberger-Treiman (GT) relation 114411 is a straightforward result of PCAC and chiral symmetry, which 
connects the ji-N pseudoscalar (Yukawa) coupling constant g„^ with the axial vector coupling of nucleon gA- Here 
its correction up to and including terms of 0(p 3 ) is obtained. In our discussion below, one can observe that the 0(p 3 ) 
correction, denoted by , is negligible compared with the 0(p 2 ) correction. 



The GT relation reads 

gAm N 



4nN - 



where Act represents the correction which can be divided into three parts 

2d 



(1+A GT ), (35) 



AG7 = _Z_^ M 2 + A (2,^ + A (3)^ (36) 
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Figure 8: (Color online) Fit up to 1.20 GeV. The fourth- and third-order fits are presented by the blue solid and green dash lines respectively. 
Results of Fit II(a)-0(p 4 ) in tableland Fit II-0(p 3 ) in table[T]are adopted for plotting. 
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loops 
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•In 

4~ Cl M 2 K 



{H (l) (m 2 N ) + {M 2 K H A {M 2 n ) - J N (M 2 )) - (m 2 H a (0) - J N (0)) + 2J?f (0)} 



(37) 



(K-2m 2 N H^\m 2 N )) + — 



m N Fl 

+ 1 -M 4 „K - \m\ (H(m 2 N ) - H m (m 2 N )) 



2(2m^ - M' K )M 2 n H (2 \m N ) - M 2 n (m 2 N - ^)£ (1) (>4 ) + M 2 A? 



Mi 



2(c 3 + c 4 )m N Ml p(3 
Fl 



2 (c 2 - c 3 - c 4 )m N M 2 „ 
H [mN) 72F& 



(38) 



The first term related to is 0(p 2 ) and generates the main contribution to Act, e.g. 1.71% for Fit l-0(p 3 ). In contrast, 
though is also 0(p 2 ), it contributes a much smaller value, ~ 0.36%. In addition, we can see from Eq. d37] i 

that its contribution is independent of the LECs. The last term in Eq. (l38l is employed to cancel the PCB terms 
generated by the terms before it. According to the naive power counting rule should be 0(p 3 ), but actually it 

possesses a chiral order higher than 0(p 4 ) and including 0(p 4 ). This can be easily observed if we reduce all the tensor 
integrals in Eq. d38l to scalar integrals and a common prefactor M* of order four will show up. It can be estimated by 



evaluating the loop integrals numerically, which leads to A 



(3) 

loops 



[-7.07ci + 1.79c 2 - 2.30(c 3 + c 4 )] x 10~ 4 . Because 



the fitted ci_4 in Fit l-0(p ) are combinations of dimension 2 and 4 LECs, we prefer to bring the corresponding Fit 
l-0(p 3 ) results in Table[T]into A! 3 ' to estimate its value, which gives a! 3) = 0.23%. 
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LEC 


Fit l(b)-0(p 4 ) 


HB^-PT 


Fit ll(b)-0(p 4 ) 


C\ 


-1.39* 


(-1.47,-1.21) 


-0.81* 


C2 


4.01* 


(3.13,3.29) 


1.46* 


Ci 


-6.61* 


(-6.14, -5.85) 


-3.10* 


C4 


3.92* 


(3.47,3.50) 


2.35* 


d\ + d2 


7.39 ±0.11 


(4.90,5.32) 


3.18 ±0.05 


dj, 


-8.04 ±0.13 


(-4.37,-3.61) 


-4.75 ± 0.04 


ds 


0.62 ±0.11 


(-1.03,-0.13) 


1.11 ±0.03 


d\\ — d\$ 


-13.90 ±0.20 


(-9.31,-8.70) 


-5.82 ± 0.09 


d\% 


-0.56* 


(-1.49,-0.84) 


-0.15 ±0.17 


eu 


3.25 ± 0.37 


(2.33,4.19) 


-9.78 ± 0.08 


e\5 


-14.50 ±0.55 


(-3.33,4.54) 


15.29 ±0.12 


e\6 


7.65 ± 0.35 


(2.74,5.69) 


-2.76 ± 0.07 


en 


8.21 ± 1.34 


(5.14,7.20) 


18.35 ±0.14 


Cl8 


-0.79 ± 1.19 


(-3.36,-1.27) 


-11.58 ±0.11 


e22 — 4e3g 


-8.19 ± 1.79 


(7.38,27.72) 


10.29 ± 0.82 


e 2 n + e 35 


-12.86 ±0.83 


(-17.35,-10.49) 


-13.12 ±0.28 


2ei<j - e 2 2 - «36 


18.18 ± 1.72 


(-25.12,-1.49) 


0.83 ± 0.55 


2e 2 \ - en 


-32.74 ±3.40 


(-7.12,-1.66) 


-25.46 ±0.48 


h A 




2.90* 


Xd.o.f 


0.03 


(0.14,0.58) 


0.11 



Table 3: LECs given by fit up to <5(p 4 ). Fit 1(b) and Fit 11(b) are performed with c, in c; are fixed at the corresponding 0(p 3 ) fit values shown in 
table[T] see the explanation in the text. Fit 1(b) is performed up to 1.13GeV, while Fit 11(b) up to 1.20GeV including explicit A(1232) contribution. 
In fit 11(b), the results correspond to c' f , d'j and e' k instead of c,, dj and et, respectively. For comparison, we provide the results from fTHl . The c,, dj 

and ek have, respectively, units of GeV~' , GeV~ 2 and GeV~ 3 , and Iia is dimensionless. The * denotes an input quantity. 



In conclusion, the correction to GT relation can be rewritten much more explicitly as 



-3.05Ji8 + 0.36 + [-7.07c! + 1.79c 2 - 2.30(c 3 + c 4 )] x 10" 2 } x 10" 2 , (39) 



The first two terms contribute about 2.07%, while the last term stands for the next order contribution which is 0.23%. 

9 (3) 

This indicates good convergence of the Act- Practically, the calculation of Act to at 0(p~) is sufficient, since A G j is 
negligible . Hence in our 0(p 4 ) fits without explicit A(1232), the parameter c/ig is fixed at the (9(p 3 )-fit value. 

3.5. pion-nucleon o~ term: cr nN 

In what follows, an explicit expression for ct„n up to 0(p 4 ) is introduced. Then fits are performed both including 
7T-N phase shift data and QCD lattice data to fix the unknown LECs related to <t k n. Finally, the fit values are used to 
predict cr„ N : cr nN -52 + 1 MeV for fit I (c) without A(1232) and cr nN = 45 ± 6 MeV for fit II (c) with the explicit 
A(1232) contribution. 

The sigma term is a quantity of great physical importance to understand the composition of the nucleon mass. It 
is defined as the matrix element of the explicit chiral symmetry breaking part of the QCD lagrangian situated between 
the nucleon states at zero momentum transfer, 



dm 



Zam N 
m„— — =< N\m u uu + m d dd\N > . (40) 
, dm„ 

q=u,d ^ 



Using the Gell-Mann-Oakes-Renner relation M 2 = Bo(m u + ma), the above equation becomes 
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Figure 9: (Color online) Convergence properties of the chiral series. The results of Fit 1(b) in Tableware adopted for plotting. The dashed 
(green), short-dashed (blue), dotted (magenta), dash-dotted(cyan), and solid (red) line represents 0(p l ), 0(p 2 ), 0(p 3 ), 0(p 4 ), and total contribution, 
respectively. 



where takes the following explicit form derived by Eq. (l2"Tt , 

m - *c lMn + ei M n - |- ^ - Ml arccos - + — In - j + ^-^ 



m N 



M$ f 3 1 M 2 3ci^M^f M 2 - 2m 2 M 2 (M 2 - 2m 2 )M n M n 
Sci - 3c3 - — C2 f In — — — \ \ ~ — In — - H - arccos 



\6n 2 F 2 n \ 4 J m 2 %n 2 Fl 1 2m 2 m 2 m 2 ■\{Am 2 ~^~M 2 2m 



(42) 



here m is the nucleon mass in the chiral limit and e\ = -2(4e22 - 8<?38 + ens + <?ii6 - 4ci?3/F 2 ). Actually, c; is the 
EOMS renormalized quantity for c\ — c\ - 2M 2 (<?22 - 4e3«). All the quantities except m on the right hand side of 
Eq. ( l42l i are substituted by the physical ones. cr nN is obtained straightforwardly, 

, , 3g 2 A M 3 n 3m 2 -Ml M n M n Mi 3 M 2 

cr wN = -\c,M 2 n + ~e x M\- * A " {— _ " arccos +M„ln-^} + — ^{8^ -3c 3 - 7 c 2 }(41n-f + 1) 

\bn l Fpn J4 m 2 - M 2 2m m \6n l Fl 4 m 2 

3c 2 M4 3c^ 2 M4 4m 2 -3M 2 12m 2 -2M 2 26m 2 M^ - 60m 4 M„ - 3M 5 n M K . 

H = — r = — =— 1 = In 1 — I : arccos — I. (43) 

647r 2 F 2 8tt 2 FI m 2 m 4m 2 - M 2 m 2 (4m 2 -M 2 )i 2m 



At 0(p 3 ) <r nN can be determined by the value of c\ . In Ref. B20I1 . through an analysis on n-N scattering partial wave 
phase shift data using the EOMS-B^PT, it predicts cr„ N =59 + 7 MeV. At 0(p 4 ) cr„ N in Eq. d43l has the unknown 
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Figure 10: (Color online) Convergence properties of the chiral series with explicit A(1232). The results of Fit 11(b) in Tableware adopted for 
plotting. The dashed (green), short-dashed (blue), dotted (magenta), dash-dotted (cyan), short-dash-dotted (yellow) and solid (red) line represents 
respectively 0(p l ), 0(p 2 ), 0(p 3 ), 0(p 4 ), A(1232) and total contribution. 



coupling constants combination e\ which does not appear in the n-N scattering amplitude. However, recently the 
lattice QCD simulations have gotten many data on the quark mass dependence of the nucleon mass, which enables us 
to fix ei as well as ci . Taking only the lattice QCD data into consideration, chiral effective field theory have been used 
to predict cr^ up to 0{p A ) II29U31I1 . In the current paper, fits are performed both including the jt-N scattering partial 
wave phase shift data and latticeQCD data. In our fits, lattice QCD data are taken from PACS-CS H], LHPC JH] , 
HSC H, QCDSF-UKQCD 51 and NPLQCD Hi collaborations. Following the strategy of Ref. Q, in order 
to minimize uncertainties of finite volume effects we only use the data with M„L > 4, and we also only choose those 
with M\ < 0.25 GeV 2 . So there are only 1 1 lattice data points which meet the requirements. They are denoted with 
stars in the tables of the Appendix A in 12911 . Note that the physical nucleon mass is included in the fits as a constraint. 
Compared with the previous fits to partial wave phase shift data, two additional fit parameters: m and e\ are included. 
The fit results are listed in Table [2] where in Fit 11(c) the leading A-exchange Born term and the partially-included 
A loop contribution (see Appendix A.2| i are considered whereas Fit 1(c) does not. The predicted nucleon mass as a 
function of pion mass is plotted in Fig. 1 1 1 1 

From Table [2] we find that most fit parameters in Fit 1(c) and fit 11(c) change little compared with Fit 1(a) and fit 
11(a) in Table|2j respectively. The prediction for cr„ N are: cr„ N = 52 + 7 MeV for Fit 1(c) and cr nN = 45 ± 6 MeV for Fit 
11(c). Our result cr„ N =52 + 7 MeV is smaller than the result obtained from the fit to the 7T-N scattering partial wave 
phase shift data up to 0(p 3 ) given in Ref. ll20Tl : cr^ = 59 + 7 MeV. We improve our determination of cr„^ in two ways. 
On one hand, the fourth order correction to cr nN is obtained. On the other hand, to our knowledge, this is the first 



20 



attempt to treat the 7r-N scattering data and lattice QCD data together using the EOMS-B^PT up to 0(p 4 ), and this 
may constrain the value of the sigma term better. The result: <t„n -4-5 + 6 MeV agrees well with the recent analysis 
on lattice QCD data using EOMS-B^PT up to 0(p 4 ) Q, which gives a nN = 43(1)(6) MeV. However, because the 
exact A-included loop graphs are not considered, cr n $ = 45 + 6 MeV can still be improved in future. 




Figure 11: The nucleon mass as a function of the pion mass. The solid line denotes the result from the Fit I(c)-Q(p ) and the dashed line is Fit 
II(c)-0(p 4 ). Sources of different lattice QCD data are: solid squares (PACS-CS) (H, open squares (LHPC) |4fj], open circles (QCDSF-UKQCD) 
Q, open triangles (HSC) O, solid triangles (NPLQCD) The solid diamond is the physical point. 



4. Conclusions 

In this paper, we performed a calculation of the pion-nucleon elastic scattering amplitude in the isospin limit within 
the framework of covariant baryon ^PT using EOMS scheme up to 0(p 4 ). The amplitude is covariant and possesses 
correct analyticity and power counting properties. The resultant description of the existing partial wave phase shift 
data from Ref. [40J] is very good for the energy in center of mass frame up to 1.13 GeV, and up to 1.20 GeV including 
the leading order A(1232) Born-term and the partially-included A(1232) loop contributions. The dimension-2, -3, 
and -4 LECs or their combinations are determined. The convergence properties of the chiral series are discussed. 
The fourth-order calculation without explicit A(1232) displays a good convergence property at 0(p 4 ) in the thresh- 
old region-the 0(p 4 ) (NNNLO) contribution is found much smaller than the LO, NLO and NNLO ones for all the 
partial waves. It is certainly an improvement to the unsatisfactory situation in the third-order calculation, discussed 



in previous literature |]23J. However, when we explicitly include the 0(p 3 ^ 2 ) Born term contribution of A(1232) in 
<5-counting 114 111 , as well as partially the A(1232) loop graphs, the convergence property is not good in the region close 
to the A resonance, which indicates that the exact 0(p 7 ^ 2 ) loop contribution may be sizable and should be considered 
carefully in the future. 

As physical applications, firstly, the correction to GT relation is discussed up to 0{p y ). The 0(p 3 ) correction is 
much smaller (about 0.2%) than the 0(p 2 ) correction (about 2%), which implies good convergence property of Act 
and confirms the applicability of EOMS-B^PT to low energy physics. Secondly, a reasonable prediction for the pion- 
nucleon <x term cr^ is obtained. We find ct„n =52 + 7 MeV from the fit without A(1232), and cr^ = 45 ± 6 MeV 
from the fit with explicit A(1232). The two values are obtained by performing fits including the pion-nucleon partial 
wave phase shift data and the lattice QCD data for m N . 
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Appendix A. A(1232) contribution 

Appendix A.l. Effective lagrangian and Leading Born-term contribution 

For our calculation here, the relevant effective lagrangian with A(1232) as explicit degree of freedom reads 

■C e ff - X.RS + -CnNA , 

£ RS = A M {iy m d a -m A r v }Ay, 

■CnNA = Tp-NT^idpKWA* +h.C, (A.l) 

where T a are the isospin-1 /2-isospin-3/2 transition matrices satisfying T\T\, - \&ab ~~ \i s abcf c - Conventions for -y*"" 



and y v can be consulted in Ref. 14111 . The leading A-exchange Born-term contribution to pion-nucleon scattering is 



0(p 3 ^ 2 ) in the so called (5-counting rule proposed by Ref. 114 ill , and the Feynman diagram is shown in Fig. IA.12l The 



\ / 
\ / 
\ / 
\ / 
\ / 

— »^ ^— » 

Figure A. 12: leading A-exchange Born-term contribution. Double solid, solid and dashed lines represent A, nucleon and pion, respectively. Crossed 
graph is not shown. 



total contribution of leading Born-term reads 



A±(s,t)=A± A (s,t)±A± A (u,t), 
B ± A (s,i)=B±(s,t)*BUu,i), 



(A.2) 



with 



B + sA (s,t) 



A~ A (s,t) 
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~2 
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3s 
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Ka(s>0 



(A.3) 



Appendix A.2. A contribution to LECs and partial inclusion of the A loop 

In order to evaluate the tree-level contribution to the LECs ci_4 when the A-resonance is integrated out, the 
leading Born-term contribution is expanded in powers of cr — s - m 2 N , Ml and f, and then compared with the 0(p 2 ) 
tree amplitude Eq. IB.4l leading to 

c a = o , £ = -d = 2c* = ^ - = l.gSGeV" 1 , 



9{m A - m N )m\ 



(A.4) 
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where we used Ha = 2.90, = 0.939 GeV, and m& = 1.232 GeV for the numerical value at the end of Eq. lA~4l We are 
not able to exactly calculate pion-nucleon loop diagrams involving A resonance. Nevertheless, this shortcoming can 
be partially remedied by substituting the 0(p 2 ) vertices in the 0(p 4 ) loop graphs shown in Fig.|4]by the contributions 
from A exchanges in Eq. iAAi . The procedure is illustrated in Fig. IA.13l Hence, the 'full' one-loop 0{p 4 ) contribution 
can be given by the diagrams in Fig. |4]with a replacement c, = c'. + c A . The effects of this replacement include only 
the 0(p 7 ^ 2 ) A-included loop graphs, while higher order graphs, like nN — > 7rA(loop) — > 7rN of 0(p ' 2 ), are beyond 
the accuracy of our calculation and therefore absent. Also, the 0(p 7 ^ 2 ) loop diagrams involving the A propagator 
contributing to the self energy of nucleon can be estimated in the same way. 

s v 
/ \ 
/ \ 

\ Q A l 

0( P 4 ) 

Figure A. 13: Matching between A-included and A-less loop graphs. The circled '2' vertex is related to c^_ 4 - Possible polynomial terms (containing 
PCB terms), denoted by a black square, would be absorbed in the c,, dj, e\ couplings as it was done in Sec. 12.51 Barring the latter, the left hand 
side A-included loop graph differs from the right hand side A-less loop graph by contributions of 0(p 5 ) and higher in the chiral expansion. 




Appendix B. Tree amplitudes 

For convenience, the two independent scalar kinematical variables, v and vg, are defined as 

s — u t- 2Ml 

v B 



Am n 4m ff 

The Feynman graphs in Fig. [2] are calculated directly and expressed in the D and B functions, while their crossed 
graphs are obtained by the following crossing relations: 

D\-v,v B ) = D + (v,v B ), D-(-v,v B ) = -D-(-v,v B ); 
B + (-v,v B ) = -B + (v,v B ), B-(-v,v B ) = B-(-v,v B ); 

A + {-v,v b ) = A + (v,v B ), A-(-v,v B ) = -A-(-v,v B ) . (B.l) 

Here the crossing relation for the A function is also displayed for the sake of completeness. Additionally, if the D and 
B functions are expressed in terms of arguments, s and t, then s should be changed to u on the right-hand side of each 
relation in Eq. IB. 1 1 Note that only the graphs (a), (d) and (e) in Fig.|2]have their corresponding crossed diagrams. 

Appendix B.l. 0{p) 
• Graph (a): 

, g 2 2m^ 

D a = ~TE^ ~ [v(v - v B ) + (m N + m 4 )v B ] , 

4b- s — mr A 

2 9 

B « = - fr-2 mN 2 (v-VB+m N + m 4 ), (B.2) 
4b A s - m L A 

where 1114 = m - 4c\M 2 - 2(8e38 + C115 + e nb )M 4 . 



9 We thank J. J. Sanz Cillero for pointing it out to us. 
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Graph (b): 



0, DT 
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2F 2 



(B.3) 



Appendix B.2. 0(p 2 ) 
• Graph (c): 



d: = - 



4c,M 2 c 2 (l6m^v 2 -f 2 ) c 3 (2M 2 - t) 



F 2 

Bt = 0, B- = 



8F 2 m 2 
2c4Mn 
~F 2 ~ ' 



F 2 



D- = , 



(B.4) 



Appendix B.3. 0(p 3 ) 
• Graphs (d)+(e): 



*1 



8aM 2 

n 

Fl 



(2d 16 - dm) 
(2di 6 - dig) 



v + 



2m N v B 



v-v B 
2mx 



v-v B 



(B.5) 



Graph (f): 



, D- f = ^ [2(di + d 2 + 2d 5 )M; - (di + d 2 )t + 2d 3 v 2 ] 



n 

4(d l4 -d l5 )vm N 



BJ = 0. 



(B.6) 



Appendix B.4. 0(p 4 ) 
• Graph (g): 

d: = 



+(e 20 + e35)M 2 v 2 + -(e 22 - 4e 3& )M 4 n j , D" = , 
B+ s = °- B ? = ^^{ei7(M 2 -^ + ei8 v 2 + i(2 e21 - e37 )M 2 



(B.7) 



Appendix C. One-loop Amplitudes 

Appendix C.l. Definitions of loop functions 

All the integrals appearing in the scattering amplitude up to 0(p 4 ) level can be generalized as 



i J (2^Aj 



#" • ■ ■ ^ 



■ A„,B\ ■ ■ ■ B„ 



(C.l) 



where A, = M 2 -(k- qi) 2 - ie and Bj = m 2 - (Pj - k) 2 - ie stems from the meson and nucleon propagator respectively. 
A standard approach to evaluate such tensor integrals has been developed by Passarino and Veltman in Ref . Q4I1 . In 
this approach, the Passarino-Veltman decomposition is firstly carried out by representing the tensor integral as a sum 
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of independent tensor structures multiplied by scalar quantities. Then the scalar quantities are further expressed by 
means of initial scalar functions of the form: 

H '""- ij (2nY'A l ---A m B l ---B„- 

In what follows, we will specify the definitions of all the loop functions and the Passarino-Veltman decomposition 
formulae, with the help of the external momenta defined as 

E„ = (P + qT = (P' + q'T, 
(T = (/" + P)f, 

*M = (q' -qf = {P-P'Y. 



1 meson: A n = H\q 



where 



• 1 nucleon: = Hq\ 



where 



2 mesons: J = H20 



where 



A A" A- -I f^LiW 
A^A^A, . j ^ m1 _^ 



Ar = rt 2) , 
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M 2 ) 2 

R+ln^-},R=- — - +y £ - 1 -ln47r, (C.3) 
fi- ) d-4 

M 2 I M 2 \ 



{A N ,A N ,A N }--J (27T)dm 2_ ( z- k )2 

K = ^a n , 



An = ^{ R+ *fj' R = jh + ^- l - ln4 «> (C4) 



Af = ^A W = ^|A W - 



m 2 \ 
32^1 ' 



1 f d d k 



) d [M 2 - k 2 ] [M 2 -(k- A) 2 ] ' 
J" Y (t) = (A"A V - g vv A 2 )J {l) (t) + A»A V J (2 \t) , 



^~\~\^ 2 -\) m 
J a \t) = i l-A* + \tJ{f) 



2t\ " 2 
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1 meson, 1 nucleon: H — H\\ 

// if /f v fl"^ = — — : - 

i J (2n) d [M 2 - k 2 ] [m 2 - (S - k) 2 ] 

fr'(s) = I?H m (s) , 



where 



H W (s) = 1 { Ak _ An + (s _ m 2 + M 2 )H( s) ] j t 

H (2 \s) = -^-y j-^ + M 2 H(s) - I(j - m 2 + M 2 )// (1) (s)J 

1 f 1 , 1 , , m M 2 +m 2 -\s' 
= 3 ~2 Aw + M 2 //(5) - - m 2 + M 2 )H m (s) + ^ 3 

ff< 3) 0) = i {-Aa, + M 2 //(s) - dH i2 \s)) , 

1 f M 2 + m 2 - -s) 
= — < -Ajv - M 2 tf(s) + 2{M 2 + s - m 2 )H w (s) ^ 3 

= -A N + M 2 H m (s) , 
H^\ S ) = l{ A ^-A^ +(s - m 2 + M 2 )H^ 2 \s)} , 

ff (5) 0) = ^- {-A;v + - m 2 + M 2 )H°\s) - 4// (4) 0)} . 



2 nucleons: 7^ = Hq2 



{Jn, 



l r d d k 
7 J 



- ; I OVui [ m 2 _ (fc _ p)2] [ m 2 _ (fc _ p,)2] ' 



where 



3 nucleons: 



= (A" A" - g vv A 2 )4\t) + A"AV®(f) + \{P»P' V + P v P'")JN(t) , 



_ i r d d k l 

H ° 3 ~ 1 J (2n) d [m 2 - (P - k) 2 ] [to 2 - (k - £) 2 ] [to - (F ~ k) 2 ] ' 
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2 mesons, 1 nucleon: 

{l,k v ,kfk v 



„ ^ If d d /t 
{H 21 ,J? 21 ,H 21 }= 7 j ^ 



[M 2 - A 2 ] [M 2 - (k - A) 2 ] [m - (P - k) 2 ] ' 



where 



^ )(t) = 5^T7{ y(f )-^ w2 ) + ( M2 -^) //2 >«} ' 

H®(t) = -^-^ {V// 21 (0 - (M 2 - 2m 2 ) H ( 2 \\t) - 1/(0 j 

= I { V// 21 (r) - (M 2 - 2m 2 )< ) (0 - i/(f)} - , 



2 [2 W \ / 21 w 2 wj 647r ; 



<« = \{\^m 2 )-\H{m 2 ) + \tH 2l {t)-H%{t) 



1 mesons, 2 nucleon: Ha 



[M 2 - A: 2 ] [m 2 - (P - A;) 2 ] [m - (P' - k) 2 ] 

A w _ ^ "A W ' 

iPf (0 = (f) + Q"Q v Hf{t) + A"A v H$\f) , 



where 



( 2 V A - 1 J „ , A _ J<f 2 rr(D, 



= I |m 2 // a (0 - M 2 H^\t) - I/ ff (0 - } , 



- - 7 {V)(m 2 ) +J H< 2) (0 



1 mesons, 2 nucleon: //^ 



{H B) J¥ B ,H B } = 7 j (2;r)d[M2 _ 



• A 2 ] [m 2 - (P - A) 2 ] [m - (E - k) 2 ] ' 

h£co = (p + 7yH$\s) + (p - syaf (j) , 

H^Jis) = g^Hfis) + (P + Ef(P + T) v Hf(s) + (P - (P - Z) v Hf{s) + 2(P"P V - H'Y v )Hf{s) 
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where 



~ 2A{ S ,m\m ^ ~ ~ ^ H{S) + K + ~ S) \ H{m2) ~ 1Mlj ^ 

+ [(s -m 2 + 2M 2 )M 2 - (m 2 - s) 2 ] H B (s)} , (C.5) 

- — \ r^r ([m 2 - M 2 + 3j1 ff(j) + [m 2 - 3m 2 - si H(m 2 ) + 2(m 2 - s)J N (M 2 ) 
2A(s,m 2 ,M 2 ) lL J L J 

+ (m 2 - s) (3m 2 - 3M 2 + s) H B (s)} , (C.6) 

~ 2-d2A( S ,Um f M ' ^ 3m2 " 3M ' + a ) Jn{m2) + K s " m2 ^ ( 2m2 " + H{m2) 

+ [m 4 - s 2 + M 4 + 2M 2 [s - m 2 )] H(s) + 2 [m 6 + 2M 6 - 2m 4 s - M 4 s + m 2 (s 2 - 3M 4 )] H B (s)} , 

\ pr- Us - m 2 ) (H(m 2 ) - H(sj) + (s- m 2 - M 2 )H m (m 2 ) -(s- m 2 + M 2 )H (l \s) 

4A(s,m 2 ,M 2 ) 1 v ; 

+4M 2 J N (M 2 ) + (0 - m 2 ) 2 - 4M 4 )h b (s) - 4(1 - d)M 2 Hf{s)\ , 

l - (2 (s + 3m 2 - 3M 2 ) J N (M 2 ) + (s - m 2 + 2M 2 ) (H(m 2 ) + H(s)) + (s + 3m 2 - M 2 )H m {m 2 ) 

4A(s,m 2 ,M 2 ) IV/ V / 

+(3s + m 2 - M 2 )H m {s) + (0 - m 2 ) 2 - AM 2 {s + 3m 2 - 2M 2 ))l B (s) - 4(1 - d)(2s + 2m 2 - M 2 )Hf{s)) 
1 



Us - m 2 ) (2J N {M 2 ) + H{m 2 ) + H(sj) + (s + 3m 2 - M 2 )H m {m 2 ) 
4A(s,m 2 ,M 2 ) 1 v ' 

-(3s + m 2 - M 2 )H m (s) + (s- m 2 )(s - m 2 - 2M 2 )H B {s) - 4(1 - d){s - m 2 )Hf{s)] 



1 mesons, 3 nucleon: 



{#13 



K) = - f- 

13 i J (2n 



H^is, t) = QfH^is, t) + (A + 2qfH^{s, t) , 



) d [M 2 - k 2 ] [m 2 -(P- k) 2 ] [m - (S - k) 2 ] [m 2 - (P' - k) 2 ] ' 

(2), 



where 



H\%,t) = - 



Ais^M 2 ) + st\ {S - m+M)HB{s) -2 {S - U)H ^ 



-{s + u) - (m 2 - M 2 ) 



Hai(t) + {s- m 2 + M 2 ) 



M 1 - -{s - u) 
2 



H™(s,t) 



^— {H B (s) - Hfo(t) + M 2 // 13 (s, t) - {Am 2 - t)H\\\s, t)} . 



After removing parts proportional to R — -jh + Je — 1 — In 4;r, the remaining scalar integrals are finite and denoted by, 
e.g. H(s), J N (f), H A (t), etc.. 

Appendix C.2. 0(p 3 ) results 

The contributions from the 0(p 3 ) loop graphs shown in Fig. [3] are displayed below, respectively. The total 0(p 3 ) 
loop contributions are given by 



A ± 

n total 



B 



total 



^ [Ag(s, tt ± A%(u, 0] + Yj A h( s > . G € {a, • • • , s) 

G // 

Y [B G 0, 0* + fi G ("> 0] + y] B„(s, t) , He {k, I, m, t, u, v} , 



(C7) 
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where the A^(m, f) and B^(u, t) (G e {a, ■ ■ ■ , i, ■ ■ ■ , s}) are obtained from the graphs (a), ■ ■ ■ , (i), (ri) ••■ ,(s) through 
crossing. 

The abbreviation F(s) in the amplitudes is defined as 



Graphs (a)+(b): 



F(s) = 2(A N - M 2 H(s)) + (s- m 2 )H {1) (s) 



Graphs (c)+(d): 



A± cd = ("2A„ + (a - m 2 ) H m (s) - 8m 2 [- (j N (M 2 ) - M 2 H B (s)) + (s - m 2 ) Hf\s)]} , 

B± cd = ■^{(&N-M 2 H(s)) + (s-m 2 )(H (1 \s)-4m 2 H i V(s)) 



i A " + ( S + 3 '" 2 ) (" { J ^ Ml ) - M 2 H B (s)) + (s - m 2 ) Hf{s))\ 



• Graph (e): 



• Graph (f): 



B} = ^{4(s-m 2 )H ( - 1 \s)+4(A N -M 2 H(s))-A n } 



Graphs (g)+(h): 



A' -A + R - -R + 



Kk = ^(s-m 2 ){-2H(s) + H«\s) + 8m 2 H«Xs)}, 

„4 

B+ g h = f^{(&N-M 2 H(m 2 ))-2A II + 8m 2 [j N (M 2 )-M 2 H B ( S )] , 

+2 (m 2 - s) - H w (s) - Am 2 (h ( b 1 \s) - H (2 \ S ))]} , 

A',, = 0, B~ gh = 0. 
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Graph (i): 

At = ^{2M 2 [H(m 2 )-H(s)] + (s-m 2 )[2H(s) + H w (s)] 

+8m 2 [(j N (t) - M 2 H A (t)) + 4m 2 H%\t) - (s - u)Hf(t) 

- (MM 2 ) - M 2 H B {s)) - M 2 (H$\s) - H i2 \s)) - (m 2 + 3s) H ( B l \s) + (a - m 2 ) H i2 \s)\ 
+32m 4 (s - m 2 )H^(s, t)) , 

Bl = {(3m 2 + s) H(s) + Am 2 H (l \m 2 ) - (m 2 + s) H m (s) 

+4m 2 [- (j N (t) - M 2 H A (t)) - 2H (2 \f) - 2 (j N (M 2 ) - M 2 H B (s)) - 2 (3m 2 + s) H^is) 
-2 (m 2 - s) H (2) (s)] + 16m 4 [(// 03 (f) - M 2 H n (s, t)) + 2 (s - m 2 ) H (2 J(s, f)]} , 

AT = --At, B7 = --Bt. 



3 ' ' 3 ' 



Graph (k): 



Graph (1): 



Graph (m): 



Graphs (n)+(o): 



Graphs (p)+(r): 



At =B + k =0, B- k = y A J<-\t) 



At = ^{4[A N -M 2 H(m 2 )]-3(M 2 -2t)[j(t)-4m 2 H^(t)]} 



Bt = 0. 



4wl V , „(3), 



AJ = —^-{s-u)H^>{t), 

2 

Bt = ~{tJ m (t)+4m z H^(t)} . 



At, = B+=0, 

A m = - 8 -^-{s-u)Hf{t), 



B- = {A„ - 4m 2 [H m (m 2 ) - (j N (t) - M 2 H A (t)) - 2H {2 \t)]} . 



At, = ^{A w -M 2 //(m 2 )} 



Bt 



2m 1 



m 2 — s 4m 



B* A ± 



2m 1 



m 2 — s 2m 
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Graph (s): 



A1 = B%=0 



• Graphs (t)+(u): 




• Graph (v): 

At = b; = o, 

Appendix C.3. 0{p A ) results 

The contributions from the 0(p 4 ) loop graphs shown in Fig. [4] are displayed below. The total 0(p 4 ) loop contri- 
butions are given by 

A fo ta , = X [Ac(s > f)± ± A c (M ' f) J + Z A±h(s ' ' G € {a,b,fl,f2,g,h,n,o,s} , 

C H 

Krai = Z[ Bc( *' f)± + B c(M)] + Z B H(^0, He{k,m,v), (C.8) 

G H 

where the Ag(w, f) and B^(u, t) (G e {a, b, f\ , f2, g,h,n,o, s}) are obtained from the graphs (a), (b), (f\), (f2), (g), (h), (n), (o), (s) 
through crossing. 
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Graphs (a)+(b): 

2c lg 2 M 2 



x ab 



F 4 



Us - m 2 )H m {s) + A N - M 2 H(s)] 



%F 4 m 2 

+2(s - m 2 - M 2 ) 



{2(5 + 3m 2 )A^ - M 2 (s - m 2 + M 2 )A n - 2(M 2 - Am 2 )A (2) 

^(s - m 2 + M 2 ) (A N - M 2 H(s) + (s - m 2 )H m (_s)) 
-(to 2 - M 2 )H m {s) + + 7to 2 - 2M 2 )H {2) (s) + (to 2 - M 2 )(s - m 2 )H (3) (s)]} 

_ (C 3 + C 4 )g (_ 1^ _ m 2 + ^2^1)^ + (M 2 _ 4^)^(2)^) + I (i _ m 2 )(j _ m 2 + ^jjQ)^ 

F 4 { 2 2 

2 

+ -^T {( s ~ '" 2 )( s + m 2 )H m (s) + (s + 3m 2 ) (A N - M 2 H(s))) 



2cig 2 mM 2 
F 4 (s — m 2 ) 



B 



ah 



2 (A;v - M 2 H(s)) + (s - m 2 )H (X \s)) 



-— ^ \-2(s + m 2 )A <2) + l -M 2 (s - to 2 + M 2 )A n - (s + 3m 2 - M 2 )^ 

2F 4 m(s - to 2 ) I 2 



H — (to - 5 + M") 

2 



-(s - to 2 + M 2 ) {2(A N - M 2 H(s)) + (s- m 2 )H m (s)) 



-2(to 2 - M 2 )H m (s) + (5s + 1 Ito 2 - AM 2 )H {2 \s) + (to 2 - M 2 )(s - m 2 )H 0) (s)\} 



(c 3 + c 4 )g 2 w 
F 4 (s — TO 2 ) 



(5 - to z + M z )H w (s) + 2(5 + 3to z - M 2 )H i2) (s) - ^(s - m 2 )(s - m 2 + M 2 )H°\s) 



c 4 g 2 m 



{4(s + to 2 )(A w - M 2 H(s)) + (s- m 2 )(3s + m 2 )H ll) (s)] 



F 4 (s — TO 2 ) 
Graphs (fl)+(f2): 

A + f = y A {(to 2 - s) [Ajv + + (s - to 2 - M 2 )H(s)] + 2m 2 (s - m 2 )H (l \s) + (to 2 - s + M 2 )H (2) (s) 

+ i(i - to 2 + M 2 ) [A N + (s- to 2 - M 2 )H {1) (s)] + m 2 (s - to 2 + M 2 )H°\s)^ 
1 _ 2ciM 2 



(C9) 



A' = -At - 



2 f F 4 

C2 



[A N + (s-m 2 - M 2 )H(s) + 2m 2 H m (s)\ 



2F 4 m 2 
1 
"4 



{s-m 2 )A {2) + -{s - to 2 - M 2 ) 



-(s - to 2 + M 2 )A„ + A^ 2) 



+-(s - to 2 - M 2 )(s - m 2 + M 2 ) [A N + (s - to 2 - M 2 )H(s) + 2m 2 H m {s)\ 

+-(s - to 2 - M 2 ) [(to 2 - M 2 ) [A N + (s- to 2 - M 2 )// (1) (s)) + 2(s - M 2 )H (2) (s) +2m 2 (m 2 - M 2 )H 0, (s)]} 



2 

C 3 fl 



B+ f 



+ |-(s - to 2 + M 2 ) \A N + (s- to 2 - M 2 )H <l) (s)] - (s - to 2 - M 2 )H (2) (s) + m 2 (s - to 2 + M 2 )H & {s) 

y A {2m [A N + A n + (s- to 2 - M 2 )H(s)] - 2m(s + m 2 )H m (s) + 4mH (2) (s) + m(s - to 2 + M 2 )H (3 \s)} 
1 \c\mM 2 



2 f f4 
+ 



2F 4 m 2 



= 1 l^^D^) + _ 2m(A (2) + A (2) } + 1 (J _ m 2 _ M 2 )C? _ m 2 + m 2 )hW 



(s) 



h,s - to 2 - M 2 ) [-6mH (2 \s) + 2to(to 2 - M 2 )// (3) (i)]| + ^2H {2 \s) + ^{s - m 2 + M 2 )// (3) (fljlO) 
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Graphs (g)+(h): 
2c lg 2 M 2 



F 4 



{(s - m 2 )H(s) + 2m 2 H m {m 2 ) - (s + m 2 )H m (s) + Am 2 M 2 (Hf(s) - H$\s)j\ 



C2g- I s-nt (2) 



2F 4 2m 2 



(A~> + Af) + -^2-0 - ™ 2 + M 2 )A n + (2M 2 - f)A A 



]-M 2 (2M 2 - t)J N (M 2 ) + M 2 (2(s - m- - M 2 ) + t)J ( ^(M 2 ) + M 2 (2M 2 - t)J®(M 2 ) 



+2(2M 2 - t)H (2) (m 2 ) + 



M 2 (m 2 +M 2 -u) + -(s - m 2 - M 2 )(2m 2 - t) 



N 

H^\m 2 ) 



+-(s - m 2 + M 2 ) [-2M 2 H m (m 2 ) + (s + m 2 - M 2 )H°\m 2 )] - (2M 2 + s - u)H w (m 2 ) 

--(m 2 - M 2 + s)(m 2 + M 2 - u)H (5 \m 2 ) + ^ - - 7 M - Us - m 2 )(s - m 2 + M 2 )H(s) 
2 4m z L 

+2(5 - 3m 2 )H (2) (s) - ((s - m 2 ) 2 + M 2 (3s - m 2 j)H m (s) - 2{m 2 - M 2 )(s + m 2 )H (3) (s)] 
-(s - m 2 - M 2 ) [M 2 (s -m 2 + M 2 )(H { ^\s) - H%\s)) - 2(2u - 2m 2 - M 2 )Hf{s) 
-M 2 (s - 5m 2 + M 2 + 2t)(Hf(s) - Hf(s)) + M 2 (s - m 2 - 3M 2 + 2t)(Hf(s) - fff (*))]} 

{(. - „r> (ifr - „? + ##> v» - fl-w) + - V ♦ « 2 - «>«"V> 

--(5 + m 2 )(s - m 2 + M 2 )H°\s) - 2m 2 [2(m 2 - s + 2M 2 - t)Hf(s) 
+M 2 (2M 2 + s - u)(H<£\s) - Hf\s)) - M 2 (2m 2 - s - u)(H^\s) - Hf(s))]} 



C4g 



F 4 

A gh = -A^(c 4 = 0). 

Ac\g 2 mM 2 



(s - m 2 ) [A„ -{s- m 2 )(H m (s) - H(s)) - Am 2 (j N (M 2 ) - M 2 H B (s)) - 4tn 2 (s - tn 2 )(H ( B l \s) - Hf(sj)) , 



j^— j J N (M') - M'H B (s) - -H w (s) + (s- m L - M')(H^'(s) - H^'(s)) 



'2F 4 n 



-(t - 2M 2 )A N - 2A™ + -(2m 2 + M 2 )(t - 2M 2 )J N (M 2 ) - 4m 2 M 2 (M 2 ) 



-M 2 (t - 2M 2 )J^\M 2 ) - ^M 2 (s - m 2 + M 2 )H(m 2 ) + ^M 2 (m 2 + M 2 - u)H m im 2 ) 
+(4m 2 + 4M 2 - t)H (2 \m 2 ) - (2m 2 + s - u)H w (m 2 ) + ( * - M ] \~ ^ 



[~H< 2 \s) 



- m 2 + M 2 )H w (s) - (m 2 - M 2 )H°\s) + hs + 3m 2 + M 2 - 2t)J N (M 2 ) 
-M 2 (s - m 2 + M 2 )H B (s) - M 2 (u - t)(H { g\s) + H%\s)) - M 2 (m 2 - M 2 )(3Hf(s) - Hf(s)) 



+(0 - m 2 ) 2 - M 4 )(H^\s) - H (2 \s)) + 2(9m 2 - s + M 2 - 2t)Hf(s) 
+{s - m 2 - M 2 )(5m 2 - s - M 2 - 2t) (h$\s) - H { ^(s)) 
-0 - m 2 - M 2 )(m 2 - s + 3M 2 - 2t) (Hf(s) - H^O))]} 



(c 3 - c A )mg 2 



F 4 



2H (2) (m 2 ) - -{s - m L + M z )H (3, (s) + -(2M L + s - u)J N (M £ ) 



1 2 -i- M 2 \tj( y U 



-M 2 [(2M 2 + s - u)H ( g\s) + (2m 2 - s - u)Hf(s)] - 2(5m 2 - s + M 2 - f)Hf(s) 

+(s - m 2 - M 2 )(2M 2 +s-u) (h ( *\s) - Hf(s)) -(s - m 2 - M 2 )(2m 2 - s - u) (lif(s) - Hf(s))) 



+Cj rf- {- 2A * + ( s - ™ 2 ) - 2H ^) + 2 (* + 3 '« 2 ) {MM 2 ) - M 2 H B (sj) 



F 4 



33 



+2(s - m 2 )(s + 3m 2 - M 2 ) (h^\s) - Hf(s))) 



B- gh = ~B + gh (c4 = 0) 



(Cll) 



Graph (k): 
At = 



(4 Cl M 2 + 3(-M 2 + 2t)J(t)\ + 2c 3 [-4M 2 A I[ + 3(2/ - M 2 ) (a„ + - M 2 ) 7(f)) 
C2 [-8m 2 A< 2 > 4- 3(2* - M 2 ) i^-t 2 J(t) + 2 (m 2 - £ ) r7«(r) - I^( f) J 



+ - 



a,: = p( w 2 + M 2 -,-£)^(f), 



4c 4 m/7 (1) (f) 



fi+ = 



• Graph (m): 



F 4 



(C12) 



3cig 2 M 2 



4 + - 



{A ff - 4m 2 // (1) (m 2 ) - 4m 2 (j N (t) - M 2 H A (t))) 



{^fV + (' - «) 2 " (3 V 2 ) - 4(2M 2 - f )H (4 V) - 2M 2 (, - ufrf^i) 



+2(s - m 2 - M 2 )(m 2 + M 2 - u) 



— (A„ - 4m 2 H m (m 2 )) - H (5) (m 2 ) + M 2 H A (t) 
4m z v 



+2f(4M 2 - t)J%\t) + 2t 2 J (2) (t) + M 2 [4(2M 2 - t)H (2) (t) + 2(s - u) 2 Hf(t) - 2f 2 i^ 4) (f)]} 
(2M 2 - t) [a„ - 4m z H a) (m 2 ) - 4m 2 (j N (t) - M 2 H A (t))} 



3£3T 
4F 4 



C4g 2 m z t(s - u)H { l\t) s-u 



F 4 



4m 



b: 



B ;, = _^|I(,- m 2-M 2 ) + i(m 2 - M + M 2 )J(^ 2 V 2 )-^ 4 V 2 )) 



BZ 



2F 4 



-J N (t) - H (V \m 2 ) + M z H A (t) - M £ H^>(t) - 2H^(t) 



A n + 4m 

+(4m 2 - t)Hf{t) - tH ( *\t)]} 



f2 w (l). 



r(2). 



(C.13) 



Graphs (n)+(o) 

,2 



Ko = -pr |(ciM 2 + 2c 4 m 2 ) [Ajv - M 2 H(m 2 )] + (c 3 + c 4 )m 2 2H (2) (m 2 ) - ^(s - m 2 - M 2 )H (3 \m 2 ) 

+-^2 [{^ ~ (ml ~ M2)2 ) { An ~ m2//(1) (« 2 )) - 8m 2 A (2) - M 2 (s -m 2 - M 2 ) [A n - M 2 (//(m 2 ) - H (l) (m 2 ))) 



B, 



2mg 



-2(s + 3m 2 - M 2 ) -(s- m 2 + 2M 2 )# (2) (w 2 ))]} 
{(2ciM 2 + c 4 (s + 3m 2 )) \A N - M 2 H(m 2 )] 



(s — m 2 )F 4 

[2(s + 3m 2 )A (2) - (s 2 - (m 2 - M 2 ) 2 ) [A N - M 2 H (U (m 2 )) + M 2 (s - m 2 - M 2 ) [A n - M 2 (//(m 2 ) - H m {m 2 ))) 



+2{2s + 2m 2 - M 2 ) (A (2) - ( s - m 2 + 2M 2 )H (2 \m 2 ))] 

+(c 3 + c 4 ) [(s + 3m 2 )H {2 \m 2 ) - m 2 (s - m 2 - M 2 )H (3 \m 2 )]) 



(C.14) 
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Graph (s): 



\ m L ) 



3g 2 Q + 3m 2 ) 
\F\s - m 2 ) 
3g 2 m(s + m 2 ) j 
F^(s — m 2 ) 2 \ 



2c x M 2 K - c 2 —A {2) - c 3 M 2 A„ 

mr 



(C.15) 



Graph (v) 

a: 

K 

Bt 



3/' 
c^(s - u)A x 



5c l M 2 A 7r - 2c 2 



3F 4 



(s - m z - M 2 )(m 2 - u + M~)A„ + A£ 



4m 2 



(21 



■ c 3 (4M 2 - f)A„ 



= 0, 



BZ = - 



Aci\inA„ 
3F 4 



(C.16) 



Appendix D. Regular parts of scalar one-loop integrals 

Following the method proposed by Refs. J6[|53l, we have derived the regular parts of scalar integrals to the order 
needed by the 0(p 4 ) calculation. All the results are listed except for the scalar integrals whose regular parts are equal 
to zero. 
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2 mesons, 1 nucleon: 
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Appendix E. Renormalization of the effective couplings 

• MS - 1 renormalized LECs: 
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Appendix F. The effect of unitarity 

The phase shift formula Eq. [33] for a perturbative amplitude automatically includes the effect of unitarity in an 
obscure way. This is illustrated in this section. 

The partial wave phase shift for an amplitude, satisfying the elastic unitarity relation exactly, can be obtained via 
Eq.[32] which reads 



,lmf' 
S f . = arctan < 



e± 



(F.I) 



Nevertheless, for a perturbative amplitude the unitarity relation may be violated. In section |231 we follow Ref. ifioll 
to define the phase shift presented by Eq. [33] Hence, starting from a perturbative amplitude, one may have several 
methods to calculate the phase shift. For instance: 



Method 1: perturbative ampl. / 



i 6 P = arctan {|^Re/ f 



Method 2: perturbative ampl. / 



p unitarization using K-Matrix method 



unitarized ampl. f K — > 5 K — arctan jf^r 



Here the K-Matrix approach is employed to unitarize the perturbative amplitude, one can also adopt other approaches 
of unitarization. f p and f k denote the perturbative amplitude and the K-Matrix unitarized amplitude, respectively. 5 P 
stands for the phase shift calculated via Eq.[33]with f p , while 5 K presents the one calculated via Eqf32]with f K . Note 
that the indices for isospin and angular momentum are suppressed hereafter. In our paper, Method 1 has been adopted 
to calculate the phase shift. Below we take the perturbative n-N scattering amplitude up to 0(p 4 ) without the A(1232) 
contribution for example to demonstrate that 5 P = 5 K . 

The chiral perturbative n-N scattering partial wave amplitude at 0{p A ) is expressed by 



f(s) = f m (s)+f w (s)+f»(s) + f m (s) 



f(2) 



(3), 



(F.2) 
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Using the K-Matrix approach, the unitarized amplitude that obeys the unitarity relation takes the following form, 



f K (s) 



/(!)(,) + f (2) (s) + Re fO)( s) + Re/W (j) 
1 - i\P\ (f (1) (s) + f {2) (s) + Re/< 3 >(*) + Re/W(*)) 



(F.3) 



where / (1) (s) and / (2) (s) stand for the contributions from the 0{p) and 0(p 2 ) tree amplitudes, respectively. The fact 
that Re/ (1) (s) = f {1) (s) and Re/ (2) (s) = f (2) (s) has been used to get Eq.|F3] The phase shift obtained through Method 
1 reads 

6 P = arctan {|^Re/ p (s)} = arctan {|^| [f (l) (s) + / (2) (s) + Re/ (3) (s) + Re/ (4) (s)]} , (F.4) 
while the one obtained through Method 2 reads 

6 K = arctan j j^^ j = arctan [\p\ [f l \s) + f (2 \s) + Re/ (3) (s) + Re/ (4, (s)]} . (F.5) 

Thus, 5 P = S K . A numerical calculation we performed also supports this observation. Hence in our work the effect 
of unitarity has already been included when performing fits to the partial wave shift data in section 13.11 The phase 
shift formula Eq. |33]for the perturbative amplitude is reasonable in the sense that it takes the effect of unitarity into 
consideration automatically. The phase shift calculated via Eq. [33] using the perturbative amplitude and the one via 
Eqf32]using the K-Matrix unitarized amplitude are the same, which is true at least for the calculation up to 0(p 4 ) in 
this paper. 

Nevertheless, the advantage of the unitarized amplitude can be shown by plotting the real part of the unitarized 
amplitude, the one of the perturbative amplitude and the unitary bound together, e.g. see Fig. lF.14l for the P33 partial 
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Figure F.14: Unitarity bound for the real part of the P33 partial wave amplitude. Left: without the A(1232) contribution; Right: with the A(1232) 
contribution. The results from table 3 are adopted for plotting. 



Appendix G. The effect of the A(1232) width 

In section [3~2l when the A(1232) is included explicitly, the phase shift is obtained through Method 1 discussed 



Appendix F 



5 P = arctan {\p\Ref(s)} = arctan {|/^ [f A ' Born (s) + f (l \s) + f 2) (s) + Re/ (3) (s) + Re/ (4) (s)]} 



(G.l) 



where f s ' tSorn (s) represents the LO Born term contribution (see Eq. IA.3) without an explicit A width in the propaga- 
tor. According to the discussion shown in Appendix F it is equivalent to calculate the phase shift via a K-Matrix 
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unitarized amplitude 
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The above equation contains an infinite resummation of Feynman diagrams shown in Fig. IG. 151 This resummation 
will generate the A(1232) width properly. 



On the other hand, one can add to the A propagator in the y A . Bo ™( s ) a n explicit width given by Ref.[41], 



r(s) = 2 



(hA_ 
\2F n 



2 s + ml - Ml 



(G.4) 



This provides us a LO Bom term A-exchange contribution with the explicit A width, which we denote it by f hY ' Born {s). 
Using Eq.[32] we have checked that the effect of f K (s) is almost the same as that of f hT > Born (s). 

Now taking +/ (2) (s) + Re/ (3) (s) + Re/ (4) (s) into consideration, one can observe that the resummation shown 

in Fig. lG. 1 Sl still exists though many new terms of different types appear. To conclude, in our work, the dominant effect 
of the A width is included through Method 1, namely through the phase shift formula Eq.|33]for the n-N perturbative 
amplitude, when performing fits to phase shift up to 1.2 GeV. 
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Figure G. 15: Resummation to generate the A(1232) width. The solid, dashed and double solid lines represent the nucleon, pion and A(1232), 
respectively. The vertical dotted lines show that the nucleon and pion in the loops are on mass shell. 
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